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ABSTRACT 


The main purpose of this manuscript is to model and 
mathematically analyze the effect of a mutualist on a population by 
its interaction with a third population. Hence predator-prey- 
mutualist and a competitor-competitor-mutualist models are 
considered. In each case conditions for equilibria are given, and 
the stability of these equilibria analyzed. As well, conditions 
giving rise to periodic oscillations are determined. Results are 


interpreted biologically, and specific examples are given. 
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INTRODUCTION 


Mathematical models in theoretical ecology, like models in other 
branches of science, are useful, because they both answer and raise 
questions. Research in this area started with the celebrated works of 
Lotka (1925) and Volterra (1931), who mostly dealt with the models of 
predation and competition. This fact, together with the literature in 
the field such as Darwin (1902) and Gause (1934) were enough to dominate 
the way of thinking of the scientists working in these areas. The 
result was an emphasis on the struggle for existence. Another side of 
nature was ignored, where there existed a variety of fascinating 
co-operative associations, even between two entirely different types of 
organisms. This pattern continues today in that most books on ecology 
pay only token attention to symbiosis, whereas full chapters are devoted 
to predator-prey, food chains and models of competition (Risch and 
Boucher, 1970). 

In its broadest sense, symbiosis is comprised of several types 
of close associations between organisms. One of the most common types 
of symbiosis, known as 'commensalism' is the relationship from which 
only one of the partners benefits, while the other is neither benefited 
nor harmed. As an example we can consider the association between 
sharks and the pilot fish (Simon, 1970). The association between them 
is commensal because the shark does not get anything from this 


association while the pilot fish are not only helped in getting food 
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from the scraps of the shark's meal but also they enjoy security 
against several predators which might otherwise prey upon them. 

The symbiotic relation between two species, in which both the 
species benefit, has been termed as 'mutualism'. Sometimes the presence 
of a third species is required before the mutualism between two species 
can be apparent. For example, in the aphid-ant system, the association 
seems to be commensal, but considering that predators may prey on aphids 
in the absence of ants, the association is clearly mutualistic. 

Mutualism is further subdivided into obligate and non-obligate 
mutualism. If the mutualistic interaction is crucial for the survival 
on one or both of the secies then the mutualism is referred as 
obligate. As an example of non-obligate mutualism we can consider the 
association between the African crocodile and the Egyptian plover, 
commonly known as the crocodile-bird. The bird feeds on various 
parasites from the crocodile's skin and the scraps of meat sticking to 
the reptile's teeth. It thereby offers a cleaning service to the 
crocodile. An example of obligate-mutualism was demonstrated by 
Limbaugh (1964) and Eibl-Eibesfeldt (Simon, 1970). It was discovered 
that the association between certain species of large fishes and 
relatively smaller fishes, known as cleaner fishes, is essential for the 
good health of the former. To know how important is cleaning for the 
large fishes, Limbaugh removed all the cleaner fishes from a few small 
coral reefs in the Bahamas. After two weeks, all except those fish 
that habitually lived in the reefs were gone and those who stayed 
developed sores and swellings. He concluded that the cleaning is indeed 
vital to the health of fish. This phenomenon is known as ‘'cleaning- 


symbiosis'. 
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Recently mutualism has received more and more attention (e.g. 
Addicott 1979, 1981; Colwell and Fuentes 1975; Halam 1980; Huchinson 
1978; May 1976; Risch and Boucher 1976; Roughgarden 1975; Vandermeer 
and Boucher 1978; Wilson 1980) as an important factor governing the 
populations of the interacting species. However, the mathematical 
theory needed to analyze such models is not yet developed. 

Kolmogorov-type models of species interaction were used by 
Rescigno and Richardson (1967) for the study of two-dimensional 
mutualism. They established conditions under which the interior 
equilibrium is asymptotically stable. Later it was shown by Albrecht 
et al (1974) that the above mentioned equilibrium is globally 
asymptotically stable in the first quadrant. The study of two- 
dimensional Lotka-Volterra models of mutualism (Vandermeer and Boucher, 
1978) has indicated that there are cases in which both populations 
survive, others where extinction is inevitable, and yet others in which 
the behaviour depends upon the initial populations. Computer simulation 
studies (Addicott, 1981) of models of two-dimensional mutualism has 
shown that these models are relatively more stable than models without 
mutualism. Basically, no general conclusions can be drawn for such 
models. Their stability behaviour depends upon the particular model 
studied. However, no model has been shown to exhibit periodic 
oscillation. 

There are many cases in which mutualism enters a system due to 
the presence of a third species. A mutualist of a prey may decrease 
the predation of its predators, or compete with its predators. A 


mutualist of a predator may increase predation on the prey, or stimulate 
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the prey to more rapid growth. A mutualist of a species may help it 
to outcompete its competitors by aiding it directly, competing with 
the competitors, or predating on the competitors. It is for these 
reasons, three-species models are important. 

This thesis is an analysis of the dynamics and stability of 
mutualism. The models treated here involve interactions among three 
species, depicting non-obligate mutualism between two of them. These 
are not just an extension of the two species mutualism to three species 
mutualism. In Chapter II, we mention some important theorems and the 
techniques which are used in the analysis of the models in the following 
two chapters. In Chapter III, we consider a model of predation in 
which there is also a mutualist to the prey. We assume here that the 
mutualist reduces the effect of predators on the prey. This could be 
accomplised in two ways: either by making prey more difficult to 
capture (e.g. Bequaert 1921; Culver and Beattie 1978; Janzen 1976) or 
by deterring the predator from feeding upon the prey (e.g. Bently 1976, 
¥577- sberger 19805) Bloomil975;; (Glynn 1976;Ross J9712 Way 1963). The 
main result of this chapter a Hopf-Bifurcation theorem, which predicts 
the appearance of small amplitude periodic solutions of the system under 
consideration, when one of the parameters passes through a certain 
critical value. We have also located the various equilibrium states 
and performed their local stability analysis. Then a special case of 
the general model has been considered in greater detail. 

In Chapter IV, we consider another model incorporating two 
competing species and a mutualist to one of the competitors. Again 


there are several ways in which the mutualist can modify the competitive 
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interaction, but the one which has been depicted in the model is by 
decreasing the effect of the other competitor. For this model also we 
enunciate a Hopf-Bifurcation theorem, predicting the appearance of 

small amplitude periodic!oscillations. Next, by taking a special case 
of the general model, we have been able to get several interesting 
results. Also, the various possible cases dealt with in Freedman (1980) 
for two-dimensional Lotka-Volterra competition models have been 
considered in the light of a third species (mutualist) and it has been 
shown that the mutualist plays a very important role, including the 


reversal of stability of the equilibrium states in some cases. 


oo - a 


1 
— ' ; ; . 
i! \ 


gee Bako: ae int bas 
ow weld eal aber ae? . Kai 


#@ soaptasqee ‘ody anki b big 
eens inisuge 2 eo Ref ik 


peony miata vit chin ‘lee curate ee ee wile: oe 
ni dud tiny 03 Fishes azipe IDV else Leste ‘dad 


dame eae #, Site cee ial Kotani ba. 4 # sae site att 


ad? giitubom: ibis rik Drocet vey o eg saiPeaeit ws 


| oo , 


BERD. ‘Biiend-f1 |, padade nag Ep he to > tala Pe: | 


rt I . 2 joe 
oe 
~~ 


1 


hy 


7 


CHAPTER II 


MATHEMATICAL PRELIMINARIES 


In this chapter we shall mention various results from analysis 
and ordinary differential equations which we shall have occasion to use 
in the following chapters. The proofs of the theorems are not given, 
but references to where they may be found are mentioned. Few 
definitions are given, but the definitions of the undefined terms and 


other relevant details can be found in the cited references. 


2.1. Characteristic equation, eigenvalues and eigenvectors of a matrix. 


Let A= Sage ie tet cee es Desa naxx pow heals Matrix. 


The equation 
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an a2 ans 


where X « ¥ (field of complex numbers), is called the characteristic 
equation of the matrix A. When expanded the above equation becomes a 


polynomial equation of degree n in i, say 


Nae ak tea fet. ."* 2: 4A + a. = 0 (Renae 


where the a's are sums of products of the ae 


The roots of the characteristic equation (2.1) are called the 


characteristic values or eigen-values of the matrix A. By the 
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fundamental theorem of algebra and its corollaries, it is established 
that there are exactly n complex numbers A, not necessarily distinct 
from one another, which satisfy the equation (2.2). 

The column vector X, = Colom ees oe is said to be an 


eigenvector associated with the eigenvalue rw if it satisfies the 


following matrix equation 


(A-A,T)X, = 6 . (253) 


where **3° 15tthe® nuliematrix: 


Zecy - Norms. 

In analysing the stability of any solution of a given system of 
differential equations, we require a measure of the distance between 
solutions at particular times; that is, of the distance between the 
vectors of solution values. Suppose that the vector function X(t) 
with projections X, (t),.--.X,(t), represents a solution vector of some 
system of first-order equations of dimension n. Several measures of 
thes sizel.cof. X(t) ) (for eee t), called norms of X, are used. It 
is denoted by ||X|| and in the simplest case it can coincide with the 


Euclidian length of the vector i.e. it is defined by the formula 


Ixtl= (3, 4) 


Two other frequently encountered norms are 


2 (2.4) 


|X|] = max |X, | Cans) 
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The norm has the following properties (shared by all) 


(CE IK || POR seomsa LT: yx 
Gray) tix|| = 0 AEP K=O 
(iii) [lox]] = fa] |[x|]], o real or complex 


|| x+¥ || <x + {|Y|| (triangle inequality) 


(2.6) 


We also require a measure of size for matrices. Let A be an nxn 


matrix oreo with real or complex elements. Then we define 


nN 


CE See 


isl j=l YJ 


2.3. Exponential function of a matrix. 


Let A be an nxn matrix. Then Be is defined by 


The series in the right hand side is convergent for all 


deduced from (2.8) that 


(i) Ie ree 


(ii) e? = I where 4 is the null matrix 
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2.4. Stability of linear equations with contant coefficients. 


Consider the general autonomous linear system 
X' = AX C20)) 


We mention the following theorem from Coppel (1965), which has been used 


many times. 


Theorem. The equation (2.10) is stable if and only if every 
characteristic root of the constant matrix A has real part not greater 
than zero, and those with zero real parts are of simple type. It is 
strongly. stable,if, and only if every characteristic root of .A. is pure 
imaginary and of simple type. It is asymptotically stable if and only if 


every characteristic root of A has negative real part. 


2.5. Liapunov function and stability theorem according to the first 
approximation. 
Consider the function V(X) 5+++5%))3 defined in the phase space 
of the variables Xp oXqee++9X)- Lope = (Xy5---5X); thens. V.0x) ers 


positive definite in a neighbourhood U of the origin if 
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If V(X) is positive definite and has continuous partial derivatives, 


then for all small enough positive c, the following property holds. 


I. V(X) < c defines an open, bounded, connected region Vo which 
contains the origin and has V(X) = c as its boundary; the 
diameter of Ve tends to zero with c; and when Cc, <c 


- boundary of Ve is contained in ve " 
if Z 


(Jordan and Smith, 1977). 
Strong Liapunov function. 


A function V(x) is said to be a strong Liapunov function for 


the system (2.10) if 


(i) V(X) and its partial derivatives are continuous; 
(ii) V(X) is positive definite in some neighbourhood U of X=0; 


(iii) V'(X) < 0 along the solutions of the system (2.10). 
Stability Theorem. (Barbashin 1970). Let us consider the systems 


T 


X" = AX +°F(X), where’ F = (F)5---+5F) C2251) 
and X' = AX (2312) 
n n 
Suppose that ‘ F*(X) ee ( Ee where o> 0 and “Ko 1s 4 
i=l 7 i=l 
positive constant. If the roots of the characteristic equation of 


(2.12) have negative real parts then the zero solution of the system 


(2.11) is asymptotically stable. 
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Domain of asymptotic stability. 

The region Rc U will be defined to be the domain of asymptotic 
stability if the Liapunov function V(X) and (-V'(X)), where V'(X) 
is the derivative of V(X) along the trajectories, have the same 


property 1: 


226. Hopf bifurcation, 

A well utilized bifurcation is the Hopf bifurcation, in which a 
fixed equilibrium state bifurcates to periodic orbits, as one of the 
parameters of the system under consideration passes through a critical 
value, termed a bifurcation value. Hopf bifurcation has found its 
applications in almost all branches of science, sometimes providing 
a useful explanation for the appearance of certain observed natural 
phenomena. 

We now state the Hopf bifurcation theorem, which is applicable to 


the types of models we are going to deal with in the following chapters. 


Theorem. (Hopf bifurcation in. R"), [Marsden and McCracken 1976]. Let 


as. be a vector field on R” for each u, Which is ck DT) OAS) 


and antl in X for each” » «(k>4)a—_ Suppose that 0) = 0 forall 
u and let A(y) = £'(0) have two distinct, complex conjugate 
u 


eigenvalues: .(:), dtu)’ “for 2 near /0, ‘such that 


(4) Re (0) = 0, F (Re Alu) |g # 9: 


(ii) remaining eigenvalues of A(y) at uw =0 have all, 


negative real parts. 


Then 


i 


ivesmmrs tn niamsh ot wa be big tie & 
YN etetie TVS bite CF « abot ve 
eniee ogtt own ic ai apts “a =e 


= 


~ : lon 


my) 


g Hotiw mi ta itoorturees gH ody ek eciacouima nt Hoek 
any to, 46 Ze taki ~iboixen oa sai kab ie i 
euisigty o gourd éoezah no bavstiaog tebe file t2 2 elt bua 
best kei dell Sebiene bd sud sulay ‘wodanemuyid'e fr 
gnlvrvetyny. wohl sites Aofngoe to aetobhye tte somite 


betuden bevigado nistise 20 gre sii? 10% vocsnbatn 
na 
J ~~ 


v — ; 
oF sltasiigqus at tikes Mo tieny, hogseouA it aqaH att othe on oie 
. 2 eee aoa ode at Ate « tgp a2 zaie wits au tobi 25. 


~ 


Sa J Loves nadsey dod bas rishi 364} Se ni ian age ; ee e 
tia) nt Aa ah dots A ins se Ny oer ibe bear be me 
lie a2. = en) Me tei eeeintue bes) “ap” tes 402 % mb Dilys b 
onmpuines xa iphion sont ‘Tib oH: ‘sued (0) 4 ms Ge ras bn co 

sda rate, 0 wea tod» im tk, ‘tounirmegia: 


70% Sahel on). Chan (Oe on te, ° nes 
ides Aral eo te) Ge 0 wewhnvodade aatareme’ (82) 
| | * eyea tone ong 


(A) There is a chs? func tan * 1, (-c =), = Ro asuch that -(s,0,u(s)) 
1s on a closed orbit of period ~ THT and radius growing 
tixe’ yi(s) for. © 7/0. “and suc that’ wt) =O. 

(B) There is a neighbourhood U of (0,0,0) in R° such that any 
closed orbit in U is one of those above. 

(C) If O satisfies a 'vague attractor' condition for the vector 


field 1e.9) CHeEMe WCS) 1 1 fOr S f0 and the periodic Orbits 


are asymptotically stable. 


2.7. Routh-Hurwitz Criterion. 

This criterion will be used to obtain conditions for the 
asymptotic stability of the equilibrium states. Consider the iit order 
polynomial equation given by (2.2). Then a formal general condition 
(the Routh-Hurwitz Criterion) can now be written, in terms of the 
coefficients Ayrsags+- 5a, which are necessary and sufficient to ensure 
that all the roots of the equation (2.2) have negative real parts. 


In our case n < 3, so that we mention the explicit Routh- 


Hurwitz Criterion;tor ni= 2° Sand. 3. 

here? a no are a OG a ee ND (25,8) 
Conditions: a, > 0; a, > 0 

n = 3 Re +a ne +a Fa. = 0 (2799) 


Conditions: “2.5% 0; a. S70: aa, - az > 0 
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the equation (2.9) has one real (negative) and a pair of pure 


imaginary roots. 


(Robert M. May, 1973). 


Zo. OUulac’s Theorem. 

The following theorem was proved by Dulac (see Andronov et al., 
Dee 205). 

Let F(x,y), G(x,y), B(x,y) have continuous second partial 
derivatives, and:let  Dyw);be Moiccig connected domain. Then if 


po) + ACB does not change sign in D and is not identically 


OX dy 
zero in any open subdomain of D, then the system x! = F(x,y), 
y' = G(x,y), has no simple closed curves in D which are unions of 


solutions of this system. 
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CHAPTER III 


PREDATOR-PREY-MUTUALIST MODEL 


S-f. General Model. 

We suppose that under the assumption of continuous birth and 
death rates the dynamics of population-growth of a predator-prey- 
mutualist system is governed by the following system of autonomous 


ordinary differential equations: 


uo = un (ux) (31a) 
KP = xg (Xx, uy = yp ee, Ww) (GSieD)) 
y' = y(-stcp(x,u)) where u,x form a mutualistic pair (SG) 


Bos < and t denotes time. 


u = the number of mutualists 
at time t 
x = the number of prey 
at time -t 
y = the number of pinaster 
at time -#t 
h(u,x) “="specific growth rate “of “the mutualist in the 
presence of the prey 
ag(x,u) = specific growth rate of the prey in the presence 
of the mutualist but in the absence of the predator 
p(x,u) = predator functional response 


Seo. @ 05 Gi>3.0 are, parameters: 
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We shall assume that 


USMY. € By 


and 


hyg,p: RXR > R, 


are continuous and sufficiently smooth functions to guarantee existence 
and uniqueness of initial value problems for (3,1) with the 

initial conditions «¢ Re and also to allow the stability analysis 

of any solution of (3.1). We require the solution to be defined on 
some interval [0,T) where 0 < T < o, 


We further make the following assumptions -- 


A(i) The mutualist (u) can grow at low densities with or without the 


prey (x). This means that for u the mutualism is non-obligate. 


Mathematically 
hes Wea aaVx. 
A(ii) The mutualist (u) cannot multiply over a certain population 
size, which depends on the population size of its partner 
ive: Bthe*prey"9(x)t* This*means*that the*mutualist**(u)» has a 
carrying capacity, which is a function of the prey population. 


This puts the following restriction on h(u,x) 
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A(iii) Multiplication of the mutualist is slowed by an increase in 


their own numbers, other populations remaining the same i.e. 
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A (u,x) su08 
This means that the mutualist exhibits density dependent growth. 


A(iv) The multiplication of the mutualist is enhanced with an increase 
in the prey population, for any population of the mutualist. 


Mathematically 
hy (u,x) a We 
This is the mutualistic effect. 


A(v) The prey can grow at low densities with or without the presence 


of mutualist, so that mutualism for the prey is also non-obligate. 
C0 sta) 0), Yur 


A(vi) The environment has a carrying capacity for the prey, which depends 


on the population size of the mutualist, 1.e. 


dk 
AK(u) a» g(K(u),u) = 0, where K: Seen asin and fae Oz 


A(vii) There may be a cost to x associating with the mutualist, 


eee 
g (x.u) <c08 


A(viii) Multiplication of the prey is slowed down by an increase in their 


own number, for a fixed population size of the mutualist. 


Mathematically 


g, (x, u) < 0. (density dependent growth) 
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A(ix) The predator functional response to prey density, which refers to 
the change in the density of prey attacked per unit of time per 
predator as the prey density changes, is assumed to be non- 


negative i.e. 


Dix,u) 220: 


Also, there cannot be any predation in absence of prey, so we 


assume that p(0,u) = 0. 


A(x) For fixed population sizes of other species, the predation is 
enhanced with the increase in number of the prey-species, that 


is 
P, (x,u) =20. 


A(xi) The mutualist cuts down the effectiveness of the predation on 


the prey. This translates into 


P,(%wW < 0. 


This is the main effect of the mutualist, incorporated in the 


model. 


These assumptions are ecologically reasonable and are as 
exemplified in nature as mentioned in Chapter I. First we establish 
that the system (3.1) is well behaved in the sense that all the 


solutions remain bounded. 


Theorem 3.1. Under the assumed mathematical restrictions on the 
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PONC CT Omcremies (idid Bp stthe#solutions u(t), x(t) and y(t) of — (3-1) 
with positive initial conditions are all positive and bounded for 


iS cage 
= 0 


Proof: Let u(t) =u X(t) = Xp and y (ty) Tek0 be the population 


Or 
Ssazes at time t= to» then we first prove that 


x Ge) < max{x,,K(0)} =M (say) 


easenit): Let x,.> K(O)e Then we claim that x(t) <x Or teat 


0 0 0° 


Suppose this is not true then 


and x'(t,) 2a 


x'(t,) = ax(t,)g(x(t,),u(t,)) - y(t) p&(t,).u(t,)) 


axe (Xa uit, ky (t p(X. .u(t, )) 


A 


axye(K(0),0) - y(t, )p(p.ult,)) (by A(vii), A(viii)) 


0 - y(t, )p(Xp-ult,)). (by A(vi)) 


SOmtidatey x(t) a0. cont rad aces On. 


Hence x(t) < Xq £O oekt > ty: 


Case. Ghigns eet XQ < K(0). Then we Claim qthat. x(t) < K(0) for all 


t > Suppose this is not true then dt, > ty 2 x(t.) = K(O) and 


> tp. 
x! (t,) 202 


We consider two subcases. 
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Subcase (1): x" (t,) 7u 0 


EEO. co). L) 


Rey) = Ox(t Jetx(t eult,))) = W(ESNP(x(tyynult.)) 
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Subcase (ii): x" (t.) = 0 


Ingthis: casesit y(t.) > 0 then we have as above 
x'(t,) < aK(0)g(K(0),0) - y(t,)P(K(0) ,u(t,)) 
=40 “contradiction. 


EE y(t.) = 0 then again we can get contradiction if g(K(0) ,u(t,)) < 


g(K(0),0), so the only case left is when 


y(t.) = 0 and g(K(0),u(t,)) = g(K(0),0). 


In this case uniqueness of solutions to initial value problems imply 
thate x<<-0 (0) lsmasolutyon. 
Thus we have proved that x(t) < K(0). Combining results of the 


above two cases we have established that 
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Case (1): Let Uy ? L(M). Then we claim that u(t) < Uy for all 


t ere Suppose this is not true then 


dt. Ss u(t) =u 


' 
> ty and u (tz) 2 MWe 


0 


Bute rrom (3.1) 


ul(t,) = u(t,)h(u(t,),x(t,)) 
= Ugh (ug ,x(tz)) 
< Ugh (LM) M) (by; A (217) vand = AGyv)) 
=O) contradiction. 

Hence u(t), <:u POs > ot 


0 OF, 


Gase..(ii9).:.. Let Ug < A(M).. nen we claim that .u(t) < L(M), for all 


Ce ot Suppose this is not true then 


0° 


Clit aoe tins u(t,) = L(M). and u'(t,) 0% 


4 0 


We consider two subcases. 


Ssubcase, (1): u' (ty) > 0, 


From: “(3215 


u'(t,) = u(t, Jh(u(t,),x(t,)) 
= L(M)h(L(M) ,x(t,)) 
< L(M)h(L(M) ,M) 


=40° ¥) contradiction, 


Subcase: (a4) u'(t,) = 0. 
Ee x(t 4) <M then again we will get contradiction because 


ul(t,) = L(M)A(L(M),x(t4)) < LO)A(L(M),M) = 0 from ACiv). 
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So the only case left is when 


u'(t,) = 0, h(L(),x(t,)) = h(L@),M) 


In this case uniqueness of solutions to initial value problems 


imply that 


Wis, (My f-rs)a solution. 


Hence: inal licases Mu(t).< LM): 


Now combining results of the above two cases we prove that 
Ut pec Nutone cre 01.0, Db. 


Finally we shall establish that y(t) remains bounded. To prove 


this we consider 


exe ty" =) claxe GG -yp(xj ule yl-step(x,u),] 
= Coxe (x;U) = sy 
= Coxe (x, uy + -cSx = scx+y) 
< AG si(cxty) 


where A = caMg(0,0) + csM. 
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Since x(t) remains bounded, this inequality proves that y(t) 


remains bounded for t > t This completes the proof of the Theorem 
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Note: Suppose that the assumption A(vii) is replaced by 
g(x,u) 2 0, 


which is equivalent to u having a direct positive effect on x, 


then again we can get boundedness results if we make the following 


restriction on the function K, 


lim K(u) = K< . 
uo 


$$ 


7h otvens) ae 


saya ewollat th aie 


aa 2 ere | “altceeay] 

{5 ; BY ? . 

de en. | : | 
0 ya, A D SB wx ae ik 
bo als 2 SE gue ho! - oftdiyr(dixo] am 

f 

(nto)? ‘ | 7 

a 23 org ee * Se Lit hee CF pea m4 


, > aa 
iH ’ yi As, 
2 = b) b 
“ity eas | . : Ve 
: fi Aor ' 
| 


gt SF METS chit git) nah & Pea), \\ 


: . ‘ it, Ue 
'S)%, Te?) aevend wiilennad: cod? , bebaaded ey Phare» (t0% 


' i, ~ : 
-g2' 37) Yo3 Gemmuads a 


i } 
movoonl alt 490 FaeTe suit ted lguigs. si fT 


9 


i} 


ve boogiqes 2% (PEr)A apie sds) vers Sa 


(D4 4 dabel a3 


ab ae aeeTie geitiade toetib Sayyid got ons Ley Lupa et aoe 
golwolle? ody aten ow 4c egigesy saanbabiuod 484 RS ow chee rindi 


2 potion why no notgisteee 
4 : : -_ we 
han - 


ee mee y ete 
one 


Ze 


woes COULLIbrium, States. 

A set of solutions of system (3.1) is the 
set of stationary solutions for which the population growth of each 
species in the community is zero. The corresponding point in the phase 
Space is known as the equilibrium point and the population sizes 
corresponding to this equilibrium point define the equilibrium state 
for the system (3.1). 

EL: (0,0,0) is always an equilibrium state. 


Promar\(21) “and “A CVI) itvisichear that 


E,: (0,K(0),0) and 


E.: (L(0),0,0) are also the equilibrium states. 
If we assume that 


= « Range p(x,0) s3e2) 
and let x be such that 


p(x,0) = (3.3) 


alu 


then 
E,: (0,x,y) is an equilibrium state, where 


ce exes) (3.4) 
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x < K(0) (3.5) 
Depending upon the number of intersections of the curves 


u = L(x) (3.6a) 


x = K(u) (3.6b) 


we will have various equilibria in the u-x plane. After Rescigno 
and Richardson (1967), we can make some further assumptions on the 
functions h and g_ so as to guarantee the existence of a unique 


equilibriym interior to the u-x plane i.e. 


Ec: (uaxje where u,Xx are such that 


and 


L(K(u) ) 
KCL) 


ul 
~ § 


Finally, we will write down conditions for there also to be an 
equilibrium interior to the first octant’ (u>0, x>0, y>0)- Any 
equilibrium of this type will be obtained by solving the following 


system of algebraic equations 


nu, x) = 0 (G3..7a) 
axe (x,u) = yp(x;u) 2-0 (2-203) 
=Sut CD (Xe) e2 a0 (S576) 


From (3.7a) we have u = L(x), so that in order to solve for x by 


the equation (3.7c), we need to assume 
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Under the assumption (3.8), the equation 
-S + cp(x,L(x)) = 0 (3.99) 


can have several solutions, giving rise to several interior equilibria. 


The y-value of these equilibria is given by 


_ axg(x,u) 
ak TCT (3.10) 


Again in order to guarantee a positive y-component, it is necessary to 


assume 
x= KCL TX) Ger) 
To have a unique interior equilibrium, we further assume 


P, (x, u) + Pb Gd =50 (S212) 


ehustunder theiassumptions | (Seciee(s29)he(seld )leandrd (34225) athere 


exists a unique interior equilibrium 


E,: (u*,x*,y*), where 


x es suchmthadL 


nou Lac) = = (34 13a) 


u’  is¥givenby 


a = L(x) (3.13b) 


and y is determined by 


noi ruta gis, (8.2) 
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* * * 
vonage ilesilae A (3.13c) 


Dix su) 


Soo, OGaDility,or Equilibria. 

The first objective in analyzing a model is to judge the 
stability of its equilibrium states. Depending upon whether the 
differential equations are assumed to apply over all conceivable 
combinations of population sizes or only in the neighbourhood of 
the equilibrium states, the stability analysis is referred to as 
global or local, respectively. The mathematical techniques have not 
yet been developed to perform the global analysis of a system like 
(3.1) in general. It has been observed that such systems can 
display a full and rich dynamical complexity such as the presence of 
strange attractors (May and Leonard, 1975). Such a complex behaviour 
is manifested even by the simple Lotka-Volterra equations for three 
competitors, which has been discussed by May and Leonard (1975). 
However, the local stability analysis can be done by computing the 


eigenvalues of the variational matrix at the equilibrium points. 


The’ variational matrix for system (3.1) 4s given by 


h(u,x)+uh, (u, x) uh, (u, x) 0 
M(u,x,y) = Jaxg (x,u)-yp,(x,u) a(g(x,u)+xg,(x,u))-yp,(x,u) — -p(x,u) 
cyp,, (x,u) cyp, (x,u) -s+cp(x,u) 
(S214) 
Now we consider the various equilibrium states separately. 


E,(0,0,0): From (3.14), we can compute the characteristic 
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equation for E,- As mentioned in Chapter II, it is given by 
h(0,0)-A 0 0 
0 ag(0,0)-r 0 = 0 
0 0 -s-A 


Me eigenvalues are “h(0;0)% a@z(050)> "=ss"“Since® “h(0,0)">" 0? “and 


g(0,0) > 0 from A(i) and A(v), the equilibrium state Ey is unstable. 


Near Ey» the u and x populations grow whereas the y-population 


declines. 


E,(0,K(0) ,0): The characteristic equation for E., is 
h(O,K(0) )-a 0 0 
ak(O)g (K(0),0) aK (0) g (KO) ,0)-A -p({K(0),0) | = 0 
0 0 -s+cp(K(0) ,0)-A 


The eigenvalues are h (0,K(0)), ak(0)g (K(0) ,0) and -s+cp(K(0),0). 
From the signs of these eigenvalues we conclude that E, 1s Stable in 
the x-direction and unstable in the u- and y-directions. Here we 


have assumed that 
S 
PULK(0) 02" mi (3.15) 


This is reasonable to assume, biologically, because in the absence of 
the mutualist and when the prey population is near its carrying 


capacity, the predators must multiply. 


E,(L(0) ,0,0): The variational matrix M(u,x,y), evaluated at 


E., assumes the form 
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L(O)h, (L(0) ,0) L(0)h, (4(0) ,0) 0 
0 ag(0,L (0) ) 0 
0 0 -S 


so that the eigenvalues are L(0)h, (10) ,0) <0 ag(0,L(0)) > 0 and 
-s < 0. This means that E. attracts in the u-direction and y-direction 


but. répels in; the x-direction: 


E,(0,x,y): This is an interior equilibrium for a predator-prey 
system, in the absence of the mutualist. We have a very rich literature 
to analyze such a system. Here we will mention a few results. The 


characteristic equation at Ey LS 


AO A 0 0 
akg (X,0)-¥p(%,0) alg (x,0)+ke, (x,0))-yp, (x,0)-2  -p(x,0)} = 0 
cyp,, (X,0) cyp, (x, 0) -A 


Expanding the determinant we get 


[h(O,%)-A] [A*-ACag (%,0) tag, (%,0)-Fp, (%,0) }+eVp (K,0)p, (K,0)] = 0. 


be i> ho and hz are roots of this equation then 
1D te h(0p,0)uacy 0 
hy + Az = ag(X,0) + axg,(x,0) - yp, (x,0) = H(x) (say) 
Aor~z = cyp(x,0)p, (x, 0) > 0 


_ axg(x,0) 


: P,, (x,0) using (324) 
p(x,0) 


where H(x) = ag(x,0) + axg, (x, 0) 
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Since the product of the eigenvalues do and re 1s positive, the real 


parts of i. and iA, have the same sign as. H(x). 


2 3 
TE eH). <= Oy E, is asymptotically stable 
Hise 0% E, is unstable (3.16b) 


in the xy plane. Freedman (1976) has given the graphical analysis of 
this case and accordingly, we can mention the following theorem for 


our system. 


Pieonenc.2. elf the assumptions’ of section gay uhold thenvatgleast one 


of the following is valid 


(re) aes is asymptotically stable in the x-y plane but unstable 
in the u-direction. 
(ii) The system (3.1) has a periodic solution surrounding Ey 


in the: x=y-plane and: lying in the strip 
MOM xy) us =0.. 80sxK(0)) yer. 


The periodic solution is stable from the outside in the 


plane but unstable in the u-direction. 


E, (u,x,0): Inthe. “u-x *planejiithe system (3.1) assumes the 
form of a Kolmogorov-type growth model. Such a system has been analyzed 
by Rescigno and Richardson (1967) and Albrecht et al. C1LOVAY inedetarl. 
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lesecp(xgu)-2) [,°-a (Gh, (U, x) +axg, (x,0)} 


touxth | (u,x) By (x,u) “hy (u,x) g (xu) iF] =a 


so that the roots i r r are given by 
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de = cp(x,u) - s 
hy + do = uh, (u, x) - axg, (x,u) <T0 
ir, = auxth, (u,x)g,(x,u)-h, (u,x)g  (x,u)t 
h (a,x) <0. g (x,u) el) ee eee 


Since ie Css ohes 0, AA, > 0, SO that in this case Ee is asymptotically 


stable in the u-x plane but the stability in the y-direction depends 


upon the sign of {p(x,u) - =}. 


Note: If we assume that the species u and x are mutualistic even in 


the absence of the predator i.e. if we assume that 


2 (xu) 2.0 (35174) 
and also make further assumptions 


uh, (u,x) - xh (x, u) Seo (Seb) 


ug (x,u) + xg (x,u) < -a < 0, (SGal7e) 


then again, we can have the same conclusion as before. Because (3.17b) 


xh (u,x) < -uh (u,x) 
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constant ratio ~ , the multiplication of each species is slowed by an 
increase in the number of individuals of both species. 
: , : uw : ; 
At Ehis “point it will be shown thats (u%x,0) is stable in the 


large for solutions initiation in the u-x plane. 


Theorem 3.3. Let h(u,x) and g(x,u) have the previously stated 
assumptions. Then (u,x,0) is asymptotically stable in the large 
for solutions initiating in the positive quadrant of the u-x 


plane. 
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and is not identically zero in R. Then by Dulac's Theorem there 
are no nontrivial periodic solutions in R. Hence the w-limit 
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(G00) 00 B: De 


; By ty Finally we compute the stability of the 
interior equilibrium. Using the relation (3.14), the variational 
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The eigenvalues of this matrix are given by the equation 
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a, = -cox*L(x*)g(x*,u*){h, (u*,x*)p, (x*,u*) -h (u*,x*)p, (xt u*)} 


using the relations, (3.13.48 ;bec) \ 
By the Routh-Hurwitz Criteria, the roots of equation (3.18) 


will have all negative real parts iff 


as > 0h, az > Oreand aja, - az > 0 C19) 


Here 
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Hence we have 


Theorem 3.4. Let the assumptions in Section A and the conditions (3.8), 
Soleo lly and “oul2jeehold= “ihenstheresexrsts a) euniquesinterior 
equilibrium for system (3.1), which is asymptotically stable (locally) 


wee a, > 0, az > 0 and aja, - a, > 0. 


Corollary 3.5. The interior equilibrium (EQ) is asymptotically stable 
if the following conditions hold 
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Proor: From (1) Wwevmget.a.>° 0, from: “<Giii)” we have 4, > 0 and 
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the-condiations., (41) “and. (iti) guarantee that - a, > 0. Hence 


le 2eaas 
by the Routh-Hurwitz criteria, the roots of the equation (3.18) have 


all negative real parts. This proves the Corollary using Theorem (3.4). 


Example: Consider the following predator prey model 


ae x(lj=x) = h2xy 


ye y Gi 2x) 


The interior equilibrium can be shown to be (% , %). For the stability 
of this equilibrium, we can apply the criterion given in Freedman (1976). 
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Hence the interior equilibrium is asymptotically stable. Now let us 


introduce the mutualist, whose effect is to cut down the predation. 


Consider the model 
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Te can be shown that, if e« and 6 are large enough e.g. ¢ = 6 = 1 
(say) then the above system does not have an interior equilibrium but 


if ce and 6 are small enough such that 
1-6 - e€6 > 0 
then there exists an interior equilibrium 


— 2t+e 1+6 (1+6) (1-6-e6) 
2-€6 2-€6 (2-e8)- 


E satisfies all the conditions of Corollary (3.5) and hence the 


interior equilibrium is asymptotically stable. 


3.4. Stability theorem according to the first approximation and the 
domain of asymptotic stability: 


We make a change of variables in system (3.1) as follows 


wo Ue = uy (S22) 
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sh aes Sade 6 (3.21c) 


where (u*,x*,y*) is the interior equilibrium of system (3.1). Now 
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Z = (u, X51) (Column vector) (322) 


the system (3.1) can be represented in the form 
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Z' = AZ + F(Z) (3.23a) 


where A = 
u*h (u*,x*) u*h, (u*,x*) 0) 
ax" gs  Cx*,u*)-y*p, (x*,u*) ae” tr txeoe(xe yu" ))-y*pi(x* ius) -p (x* ,u*) 


* e ak 
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and F(Z) is a column vector given by 
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The linear approximation of system (3.1) near the equilibrium state 
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Z' = AZ (524) 


Now we state the following theorem. 
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Theorem 5.0. Let the assumptions in Section A and conditions (3.8), 


(Seo )y> (S-91)- and “(3.12)) hold, and let a,24,,a, be as in (3.8). 


Also let 
(a) a1 = (le a > 0, “and aja, - az > 0 
ae Ss WECZ) Nt 
So ola 


for any t, > 0 is an asymptotically stable 


Piverer a(t) = 0) ts oot 0 


0 
Soiution of (3.23) 


Broor-> We outline the proof, which is given in Jordan and =smith (1977) 
because from that we would be able to say something about the domain 
of asymptotic stability. 

Due to condition (i), A is a stable matrix, so that there exist 
constants  c > Ose ye< Oeesucnstnat 
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We can construct a strong Liapunov function 
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ViQZa ie tke (3.26A) 
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The condition (3.25) ensures the convergence of the above integral. 


V(Z) is positive definite, which becomes obvious when we rewrite it in 
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Here the integrand is simply the sum of certain squares. 
the time derivative of V(Z) along the solutions of (3.2 
the relation (3.26A). 
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Now consider the product (e ). We can write 
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which gives 


-I = alx + KA. 


Using this relation, the equation (3.27) becomes 
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Now we compute 
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This proves that Z(t) = 0 is an asymptotically stable solution of 


C3525). 


Domain of asymptotic stability: 

From the above proof it is apparent that it has not only been 
proved that the zero solution is asymptotically stable, but that all 
solutions starting in certain neighbourhoods of the origin are 
asymptotically stable. Such a neighbourhood determines the domain of 
asymptotic stability, with respect to the particular Liapunov function 
considered, so that this does not necessarily determine the largest 
domain. Now as mentioned in Chapter II, the domain of asymptotic 
Stability: LR(Vie)) will be taken to be the region in’which V(Z) ‘and 


(-V'(Z)) have the property I. So that we can say that 
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Soja POriodic Solutions. 

In this section we obtain the conditions on the parameters to 
guarantee the existence of small amplitude periodic solutions for the 
general system (3.1). ° As mentioned earlier in Section C, in stepping 
up to three dimensional models, numerous problems are encountered. 
Theorems like the Poincare-Bendixon Theorem do not hold for dimensions 
greater than two. In such a situation, the Hopf Bifurcation Theorem 
1s one of the important direct tools, which helps us to establish the 
existence of the periodic solutions to non-linear systems of order 
greater than or equal to two. Experimentalists have shown that sometimes 
oscillations occur as a result of changing some condition or parameter 
in the environment (see Bunning 1973), so that theorems of the type 
mentioned in this section could be of great help in providing the 


mathematical description of these phenomena. 


Let (u*,x*,y*) be the interior equilibrium state for the 
system (3.1). The characteristic roots of the variational matrix 
evaluated at (u*,x*,y*) are given by equation (3.18). Rewriting - 


this equation we have 


a ae + ash + a, = 0 (3.29a) 
where 
ceo eax ee oa) YD Ge ou (35296) 
a, = [ag(x*,u*) (cx*p, (x*,u*) +u*h, (u*,x*) ) 


+ou*x*(h, (u*,x*)g, (x*,u*) hy (u*,x*) g (x*,u*)) 


Psy ACW (Hemel petats a yume paces) (3.29c) 
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and 
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Now we state the following Hopf Bifurcation Theorem for the general 


model (3.1). 


Theorem’ 3.7... Let. (u~,x* vy") 9 be the interioriequiiebrium statesat, the 


system (3.1). Further, let the following conditions hold 
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appear, which bifurcate from the equilibrium state (u*,x*,y*). Here 


@),45,a, are given bY wi(3720) wand by »b,, DYF oS 230). 


Proof: First we show that when a = Oa» the characteristic equation 
(3.29a) has one real and two pure imaginary roots. For this we compute 
(a, a,-a, and express it as a polynomial in a. Using the relations 


Cor) tS. 29C)} ands. Co. 29d) are also (>. 15c),, 1t can be shown that 
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Nowousing (3213a)\, (6.30a)) Gand (6-30b), the above expression can be 


rewritten as 


Hence when a =a 


= ap (byagtb) = 0 (3.31) 


(a,a,-a i%0 


3)a=a, 
The condition (i) together with the relation (3.31) guarantees that 


when a = Cty» the characteristic equation (3.29a) has one real 


(negative) and two pure imaginary roots. 
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Next we show that the eigenvalues cross the imaginary axis with 
non-zero speed at a = Og - We can assume that for the values of a 


near a » the characteristic roots are of the form 


and are real numbers. In view of this the 


where hyo ho 3 


characteristic equation will be 


24s 
(A-AZ)LO-A,) tpl = 0 


or 
3 D 2 MZ Zee . 
ie (Az+2A,)A + (Ap tA,FZAVAZIA - (A, tA,)A2 = 0 (3.32) 
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Now let us differentiate this relation with respect to the parameter a 


throughout. We get 


da da Ghee \n eda di 
2 1 nea ee 3 1 2 
da ay * 2%, ) fi ii a4 da yi da : da plea 
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~ 20, °2(a,+2X,) (eee 2 sane 
(3730) 
Now consider the equation (3.36) for o = Og - Since we know that 
at a= Oy» two roots are pure imaginary, 
(Biles & = 0, (335375) 
0 
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0 0 
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ag 0 
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Hence atl is We nd # 0 (3.38) 
GOT ees 2 a“+a 
OG) a A=. 


This shows that the eigenvalues cross the imaginary axis, transversally, 


i.e. with non-zero speed. 


Now the application of the Hopf Bifurcation Theorem, given in 


Chapter II, proves the theorem. 
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2To.h ie SpeClalmGase: 
In this section we analyze a special case of the general model 
(3.1), incorporating all its important features, as mentioned in 


Section A. We consider the following model 


Bie raf aa eh 
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x' =ax(1-s)- 
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whererthe parameters, a58.7,%,l.,Kk,m,¢,s sare all positive. This 


0? 
particular model refers to the case, where in the absence of the 
predator (y), the association between u and x is not mutualistic 
but commensal. The mutualism occurs when we introduce the predator into 
the system. 


The interior equilibrium is obtained by solving the following 


system of algebraic equations 
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Putting this value of u in the first equation; we get 


s(1+mL,) 
5 cB8-2ms ” 


and then finally from the middle equation we solve for y: 


1+mu 
8 


Doh) 


We denote the interior equilibrium by E*(u*,x*,y*), where 


us = Ly 1 = r 
x* = = (3.40) 
y* 2S eh 
where 
cB (1+mL) ) 
L = Pee he (Sorel) 


To make the equilibrium E* feasible, (i.e. positive), we assume 


cB > &ms (42a) 
and jxaee’ giv (3.42b) 
BK 


The latter condition requires that x* < K, which is the usual assumption 
in two-dimensional predator-prey systems. 


Other equilibrium states are: 
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EL: (0';,.0,50)) 
E,! (055630) 
ER: (L,,0,0) 


ee (pee 
Ee (0, eaeia bate 


Eo: (Lot2K,K,0) 


We will concentrate on the interior equilibrium E* only. The 


variational matrix for system (3.39) is 
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The eigenvalues of this matrix are given by the roots of the following 


cubic 


He + au wea + ae =U (S44 a) 


where 


, ieee 
Th aes 


a (0,8. Ase) “3 
SE ahah rita oa 
; -) 
sai. “ino “2 «muiddci ies retrermu tay NG» nee 


aE (Chea, nierexe 103 sinsoe 
ee ee 


j P ae 
ine aD 
es ngit) Ay 
‘ume sim f A 
Pits ba 
{owen *" Pern 
7 


26g oW (00.7): prbew . (%, *%, eee at 
_ i 7 ; ‘ a 


4 # 


. . s : hs 2 . Gane Hest = : : 
(28.7) - one | , ee a) cae = un aE . 


a ae a) n sy 7 : c 
‘ y, 3 bd ~ . Cy 4 a é ! 
ied 0 ‘ete = 47 Soe oT) | 
7 - ; ¥ J 7 2 J - - =" : ? : 7 
“ Py _ e- * 


48 


B asa 

A A bea a 
7 2 asyA 

a, (1 - Sa ——)as + aK. (344¢c) 
eae 

ag = mr: i - $*) (c8- gms ) (3.44d) 


Clearly ana 0 and from? +(3:.42) a > 0. Hence the roots of 


equation (3.44a) will have all negative real parts iff 
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The above discussion, proves the following result. 


Theorem 3.8; Let the. conditions (3.42a) and (3:42b) shold.) Then: 
there exists an interior equilibrium of the system (3.39), which is 


asymptotically stable (locally) if the parameters al@.y,0,.u,,4;m,c,s 
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Proof: The conditions ms < c8, sid < c8K guarantee the existence of 


the sinterior equilibrium !r*, ras given by (5.40). From: (i)) we have 


&my gms ae & : y 2m a3 
cB le aa Pandey Ss this implies that Oe rare 202 Prom. ((45) 
we have a- > O° * Now a-S PO) Pa -S 7) Mande se. 7%. 


imply that the condition (3.45) is fulfilled. Hence the characteristic 
roots of the variational matrix for the system (3.39) at E* have all 


negative real parts. This proves the corollary. 


Fie In this section we estimate an existence region of asymptotic 
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Here u*,x*,y* are known in terms of the parameters and are given by 


(3.40). Now let us introduce the column vector 
aly 
Z= (u, x) 7) 
Using this the above system can be rewritten in the form 
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where A is a 3x3 matrix given by 
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The linear approximation of the system (3.39) near E* is 
Za 2 hZ (3.49) 


BS’ discussedteariier, A). is ‘a stable) matrix, if the}condition 4{3.45) 


is satisfied. 


Liapunov function: If Z(t) = 0 is asymptotically stable then according 
to the standard ODE thepry there exists ’a Liapunov 


function for (3.49). -The Liapunov function can be taken in the form 
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Equation (3.50), then becomes 
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From the third of equations (3.53) 
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The other elements of the matrix B 1. by Py 29) > 22532 are 
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Using the Cramer's rule, we can solve the above system if the determinant 
of the coefficient matrix is non-zero. We show that the determinant 


of the coefficient matrix, denoted by A , is non-zero by evaluating 
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Since we have assumed that the real parts of the eigenvalues of all 
the roots.of the matrix A are negative, the condition (3.45) is 


satisfied and hence 
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Thus we can solve for By 2 2by 32 59b 35 from, (3255). using. Cramer's 
Rule 
-1 2a.) 2a.) 0 0 
Bt ae a Ont cae ee 
Lata fs a ais omy sa ae ES. 
‘ SB 333 PTT dine 231 
A 2 “ Te eta es? 
2a5) -1 2a.) 0 0 
Voibeecale = oh? sBzane aan 
Say ea aes 220 eee oD) 
C : 2 a0 hae al 
F ; ca, ’ ; . 


(me. é% 
; ; 7 7 
ile Yo eeutevnoyia ot Yo arteq fast sas tant): Goins. vt 
7" a = ae oe 
el, (88,8) aelr2betes. egg eviesgen es6 & dese othe 1 S500" 
' soined bite . it : 
(8.2) : O<5 


S 


ay j ' ; 7 
2* temes?. abet, (SCcee) ~~ mots ar © cetony drgpden et ao wise mAs shal 


as 


” o» fee Leds 
pO hea 0 Sh 
(58c.¢) | ' _ <2" 
. | tc? . eee pe 
. i v » on m 
|) fe i ‘6h 
hte » ee a 6 


(Gte.£) | a ie 


% 
ian 
se ay Ct 


56 


2ajy 2a.) - 1 0 Q 
0 Bar st PBs lis =e 2ay.y, #0 
Pas F = dv2eeanl® 322 ; as Joy |S (3.58c) 
p a2 F =i Eriol 
: == £23 enans2 
2ajy 2a.) 2az, - 1 0 
0 Dae) 0 ae 0 
Cn = 5 Os bre Macy, a (Oar 
0 a a ~ maz a 
ey 11 25 31 
: u 210 ; a 252 
2a) 2a,, 2a.4 0 -l 
0 eae 0 2atee eee ae 
aa) : ieee ieee oso” Bo ; os aS) 
0 a3 hi tae ) =i! 
y : Sipe ok ; a 


Thus? the relations ~(6.52)5 (754) ($256), (3.584,b.¢,d-6). determine 


the matrix B completely in terms of the parameters of the system (3.39). 
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The function V(Z), then can be expressed as 


ED by 8) J Boa ‘ bss¥ ADE TRS Benoa eae Mace) 
We prove the following result (E.A. Barbashin 1970) 
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Proof: Computing the derivative of V(Z) along the solutions of 


(3.46) we get (Barbashin 1970) 
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Since |V'{Z) <0 in WZ <k, 4{-V'(Z)} has the property Tin WZWe<k. 
The interior of the™ball iZll < k- in-which _V(@Z) has the property I, 
will then determine the region of asymptotic stability. This domain of 
asymptotic stability depends upon the particular Liapunov function 
chosen and hence does not determine the largest such domain. 


This proves the theorem. 


Theorem 3.11. (A Bifurcation Theorem). Let the parameters a,8,y,2,L9, 


K,m,c,s, all*positive be such that 
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then as the value of a (bifurcation parameter) passes through Oo » 
we get 'small amplitude' periodic solutions of (3.39), which bifurcate 


from the equilibrium state E*, given by (3.40): 
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Proof: The condition (i) and (ii) guarantee the existence of the 
interior equilibrium E*, given by (3.40). We will use the Hopf 
bifurcation theorem for vector fields to establish the appearance of the 
periodic solutions as one of the parameters a passes through a critical 
value Oy - 

First let us show that when a = Oy» the spectrum of eigenvalues 


at E* has one real (negative) and a pair of pure imaginary roots. The 


characteristic equation at E* is given by (3.44a)* ‘i.e. 


u +a, = 0 (3.60) 
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or 
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Now as in Section 3.5, we can prove that 
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Now applying the Hopf bifurcation theorem and (3.63), (3.64), the theorem 


is proved. 


The stability analysis of the periodic orbits can be done using 
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messy that we have deferred them to the appendix. 

As a numerical example to illustrate the above, we take c = &% = 
Weees = lL = 1, B= 25°) = 4, 1k = 8. “With these values, conditions 
(i), (11) and (iii) are satisfied and the value of dy comes out to 


be 8. Hence if 8 is in the range of values of a, the equilibrium 


state (3,2,5 a) bifurcates into periodic orbits. 


oeec, tne. Case or No iInterion Equilibrium. 


Here) We consider what happens when there is/no interior 
equilibrium. Before we discuss this we need to know about the local 
stability behavior of other equilibrium states, i.e. ELE, Es.Ey and 


E. of the system (3.39). 
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S 4 ; ; 
Cor ee unstable in the y-direction 

Ss : : ; 
CorS Co stable in the y-direction. 


The latter case, in which E., is stable in the y-direction, shows that 


even if there is an abundance of prey the predator population declines. 


This could happen if the predators are unable to capture prey or there 
is some other agency which deters the predator from feeding upon the 
prey. It can be shown that in this case the predator population goes 


forextinction.. .From | (3.439) 
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so that y(t) >0 as t+. Extinction, as modelled by our system 
(3.39), cannot occur in a finite time. We will show later that in this 


particular case, the system does not permit an interior equilibrium. 
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E,(L,,0,0): Thevvertational matrix 2c. 6 is given by 
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The eigenvalues of M(E.) are -7y,0a,-s. Thus E. actracts in thes u; 
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asymptotically stable in the x,y plane. . We also note here that 
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The eigenvalues of M(E.) are 


sas aur -S iy. SCE Naess 
1+m(L)+2K) 


so that E. isestablesingthes u,x. directions but the«stability, in 


the y-direction is determined by the sign of the third eigenvalue 


1+mL, 
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cB.o> sas + ) => unstable in the y-direction 


1+mL 
CB = stmt + => stable in the y-direction. 


Now we consider the circumstances, in which there is no interior 
equilibrium. From (3.40), we find that there are two cases in which 


tats, can occur. 
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Now we show that in both the above cases the equilibrium state Ee is 
asymptotically stable in all the three directions. For this we need to 


show that 


cBK 
ui Tem(L5*2K) ow 


Consider the case I 


oak -s (1+mL,)-misk+cBK 


7 em (Lara) T+mL )+m2K 


-s(1+mL,)+K(c8-mgs) 


1+mL,+m&K 
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In the case II 


Ag Oi Sy eset Be sy an 
K K cB-2ms K \c8-2ms 


multiplying both the sides of this inequality by K(c8-%ms) > 0, we get 


(c8-2%ms)K < s(1+mL) ) 


This again implies that 


cBK 
¥ T#m(L*2K) : 


Thus we have established that if there is no interior equilibrium then 


Ee is asymptotically stable (locally). 


Theorem 3.11. If cg < z then the equilibrium state E, (L)+2K,K,0) 
is globally asymptotically stable in R° : 

Proo£s.+ «Let (wet p2cGt isnt) ) representa solution ‘to, (3.39): with 
initial conditions u(t) = Up; X(to) = Xq> y (to) TiN ge defined for 


aieist > to: Then 


moO lee yee soc: 


is a semiorbit of (3.39). As established earlier in Eo. we can prove 


that in this case 


lim y(t) = 0 


to 


so that L(c’) lies in the u,x plane. The theorem then follows from 


Theorem 3,3. 
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Observations: In the’ absence of the mutuabist (or (m = 0), the 


predator prey system is 


x! ax(1 - 2) - Bxy 


yo = sy,(-S+c8x) 


Using the notations given in Freedman (1980), the sign of 


x*g(x*)p, (x*) 
HICKS Een" BOK") tye RS) de rea 


decides the stability of the equilibrium position 


as 
cBK 


Stable. Now we want to show that when the mutualist is introduced into 


in anaes. case. H(xtj}i-<s- < 0, so that the interior equilibrium is 
this system, existence of the interior equilibrium depends on the 
parameter lomicinisuchcatiway, thatithere existseaeritical value, m* of 
m given by 


c8K-s 


SL)+%SK 


m* = 


such that if m > m* then there is no interior equilibrium. It can 

be shown that m > m* = y* = FAG -3) <0. for thee values sof  me< m*, 
we always have an interior equilibrium given by (3.40). We also 

observe that if m is small then it has a positive effect on the whole 


community in the sense that equilibrium population level of all the 


species is enhanced compared to the case when there is no mutualist. 
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Calculating the partial derivatives of u*,x*,y* with respect to the 


parameter m, we find that 


gu* _ a) GS 
om cg ‘om 
gO een 
om ch ~om 
oY Phasts _ 28h4 (9A. 
om Ss B (t ak 3m 
3d _ cB (cBL)+£s) 
where cri a aie 0 
(cg-2ms) 


These relations show that u*,x* populations increase with the increase 


of m but the change in y* depends upon the sign of the factor 


2sr 


(1 - ay: Itican ber shown thatwif 
CBK=2S. py 2sX 
m < ZsL)+isk = (m, say) then 1 - aK 2 0; 


Thus we find that as long as m<m, (< m*), the interaction with the 
mutualist is helpful to each species in the community. A similar 
result has been obtained by Thomas G. Hallan (1980), while studying the 
effect. of cooperation on competitive systems. But if m, <m< m* then 
we observe that the equilibrium population level of the species u and 
X goes up but that of y declines. This means that there are 
many mutualists protecting the prey from being attacked by 
predators and thereby causing non-availability of food to predators in 
spite of the abundance of the prey. 

In the case c8& is very small, we have already shown that 
predation is so poor that the predators cannot survive and the interior 


equilibrium in the u,x plane is a global attractor in R? 
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As regards the stability of the interior equilibrium (Ue Ne 
we find that it continues to be stable as in two dimensions (1-4ece 1 
the absence of the mutualist) if s 2 ys (Gorollary 3.9)= Thus) the 
relative magnitude of these parameters is important. We also find that 
the growth rate of the prey species (i.e. a) plays a very important 
role. When y > s, then there exists a critical value bg of as such 
that if a goes below a » then there is a change in the stability of 
che equilibrium state. o = oh becomes the bifurcation point and we 


fing that. close to Oy » the system exhibits periodic oscillations. 


3.8. Summary. 

In this chapter a predator-prey-mutualist system has been modelled 
and mathematically analyzed. Conditions for equilibria were given, and 
the stability of these equilibria determined. Conditions were also 
given for the existence of three-dimensional periodic solutions. A 
specific example was discussed. 

It was found that by adding a mutualist to the system, the prey 
equilibrium value is increased. This has the effect in the case of a 
stable interaction of increasing the effective carrying capacity for 
the prey. Further the carrying capacity for the mutualist is also 
increased. However, as expected, all populations remain bounded. 

Depending on the parameters, adding a mutualist to the system 
could be either stabilizing or destabilizing, and therefore limit 
cycles could appear, where they were not before, or disappear. 

Finally, adding a mutualist to a predator-prey system could 
cause the predator to go extinct, in which case the prey and mutualist 


population numbers approach equilibrium values. 
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CHAPTER IV 


COMPETITOR - COMPETITOR - MUTUALIST MODEL 


Study of models of competition was started as early as 1920's 
by Lotka (1925) and Volterra (1931). The models considered were just 


the extension of the single species logistic equation: 
X) = y,Xx (ee Seen ee )] 
}. bel K le Slee 92 
eee x eos Soe Pe 
2) ie heals ate ic 


where are «the intrinisic¢ jgrowth rate, «carryingscapacity .and 


ee iia 
the competition coefficient (which measures the inhibitory effect of 

xX, upon x); respectively, of the species Xy> and ¥ 2K, 205) are 
Similar parameters for the species X, + The stability properties of 
this model are well known. In particular, when the two species use the 
resources in the identical fashion (whence Oy 5 = 5) = iteand K = K)3 
2t*1s"known. that both'ispecies ‘cannot: persist. This led to the 
'competitive exclusion principle', that species which make their livings 
in identical ways cannot stably coexist in a stable environment (May, 
1976) is confirmed in the study of simple laboratory systems. Gause 
(1934 )and others, e.g) Levins) (1974), May (1973,41976), Prelou.(1977), 
Roughgarden (1975) have discussed this problem. 

Several more realistic models of two competing species have been 


analyzed by Hutchinson (1947), Cunningham (1955), Utz and Waltman (1963), 


Rescigno and Richardson (1967). 
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These models can be improved by introducing a third species in 
the system, which for example might help the weaker species to coexist 
with the other species. Hallam (1980) has discussed the effects of 
cooperation on competitive systems, modelled on the basis of Lotka- 
Volterra kinetics. He considers a system composed of a competitive 


subcommunity and two cooperative subcommunities and concludes that 


mutualism can be beneficial for all species in the community. He also 


shows that a mutualistic species can drive a stable competitive 


system to extinction. 


In this chapter we shall analyze a model incorporating two 


competing species, Xs Xy and. a third spécies wu, which acts vas a 
mutualvst to the’ species xy - We shall assume that there is no direct 
interaction between u and x,. 


2 


wae pGeneral Model . 
We shall suppose that the dynamics of population-growth of 
a competitor-competitor-mutualist community is represented by the 


following system of equations: 


u' = uh(u,x,) 


x4 rz x5 [85 (x5) -4,,(X) »X5) ] 


where u, x, form a mutualistic pair. The hypotheses implicit in 


i 
equations (4.1) are that the rate of increase or decrease of the 
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populations does not depend on time and that the populations are so 
large as to be measurable with real numbers and not subject to random 
fluctuations. 


We shall assume that 


and h, 8) 85545: Re eae: R° 2 aks 


are continuous and sufficiently smooth functions to guarantee existence 
and the uniqueness of solutions to initial value problems for (4.1) 
With Initial conditions in hee and also to allow the stability 
anabysise0t, any. solution Of, (4.1). -We require the solutions, to be 


defined on some interval [0,T) where 0 < T < ~, 


We further make the following assumptions: 


(i) h(u,x) satisfies conditions A(i)-A(iv) mentioned in 
Chaptersi1l, 
1 1) g, (x,u) satisfies conditions A(v)-A(viii) mentioned 
inscGhaptersLL I: 


(jae, The especies xX, can grow at low densities i.e. 


g,(0) > 0. 


(iv) The environment has a carrying capacity for the species 
3? i.e. 
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(v) Multiplication of the species x, is slowed down by an 


increase in their own number. Mathematically 


2 ono 


1.€. the species X, has density dependent growth. 


(vi) In the absence of the species x there is no competition 


2 


€O species x so that its growth is not inhibited. Mathematically 


1, 


q, (x, ,0,u) =, () *tor all X,.u > 0. 


(vil) For a fixed population of x,, an increase in the popula- 


‘V2 


tion of the species xX, inhibits the growth rate of Xy> l.e. 


ee Oleg > 0 


This is the competitive effect of xX, on X,. 
(Vill) For a fixed population of the species xy and X5, an 
increase in the population of the mutualist species u, reduces the 


competitive effect of xX, on XxX). Thus mathematically 
CRS Boe 


This is the main mutualistic effect of u on Xi 
(ix) Other variables remaining fixed, an increase in X1> 


results in more competition, i.e. 


q > 0 
1x, 
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(x) In the absence of the species xy there is no competition 


EO eSpecies X55 ne. 


q5(0,x.) = 0 for Xo > 5 Oe 


(xi) For a fixed population of the species X5, an increase 
in the population of the species x) Slows down the growth rate of 


% i.e. 


Diz 


nae? ae 


This is the main effect of competition of X, On X,. 


Carr) Ix, S105 


(xiii) The functions q, (X1 9X5 0¥) and q5(X, >X5) are non-negative. 


4.2. Boundedness of Solutions. 
In this section we shall show that conditions assumed in the 
Section 1 ensure the boundedness of solutions of the system (4.1). 
Let u(ty) = Up» x, (to) = Xj and xX, (tp) = X5Q ? 0 and 
represent the initial population of the species at any time to: 


Then we first prove that 


x, (t) < max{x K, (0)} 


10’ 


Caser/(Gs).;) “Let X19? K, (0). Then we claim that x, (t) < X19 for 


t >t). If this is not true then dt, > ty such that x,(t,) = x, 


and x; (t1) > 0. But from (4.1) 
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Ky (ty) = ox, (t,) [e, x (t,) .u(t,))-a, Gx, (t,),x,(t,) u(t,)0] 


= OX) gl8y yg -4 (ty) )-4; yg, (t,) u(t, J] 


[A 


AX o lS, (X1 > 0) my (X, 9X5 (ty) wee, 0] 
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QO” from *(xXiit)} and the faet. that g, (x Op aces) 


LO: 


contradiction. 


Hence x, (t) < Xj: 


Case (ii): Let X10 i K, (0). 


Then we claim that x, (t) < XK, (0) FOr all teat Tienot. 


0° 
= ' 
then dt. > 0 such that x, (t,) K, (0) and xj (ts) > oe 
We consider two subcases. 
Subcase (i): nee won 


From (45 1) 


x1 (t,) = ak, (0) [g, (Ky (0) u(t,))-a, (K, (0) .x, (ty) u(t.) )] 


[A 


ak, (0)g,(K,(0),0) (by (11) and (xii1)) 


O; contradiction: 


Subcase (il): x; (t5) = 0 
In this case if g, (K, (0) ,u(t,)) < g, (K, (0) ,0) then again we 
can get contradiction. If g, (K, (0) ,u(t,)) = g, (K, (0) ,0) and 


q, (K, (0), x, (t,) ,u(t,)) > 0, then 
ey = -aK, (0)q, (K, (0) .x,(t2) u(t.) <50), contradiction. 


So the only case left is when 
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g, (K, (0) ,u(t,)) = g,(K,(0),0) and q, (kK, (0),x,(t,) ,u(t,)) = 0. 


In this case uniqueness of solutions to initial value problems imply 


that xy = K, (0) is a solution. 


Thus we have proved that x, (t) < K, (0). Combining results of 


the above two cases we establish that 


x, (t) < max{x,),K, (0)} = Mo itsay)i 


In a similar way, we can show that 
X(t) < max{X5.,K, } =N (say) 
Finally we prove that 


u(t) < maxtu,,L(M) } 


Let Uy > L(M). Then we claim that u(t) uy, fOreal leet pat 


Suppose this is not true then 


0° 


= t 
dt. ae Oar u(tz) =u and u (tz) PO: 


0 
From, (4.1) 
ba ike 


= u(t,)h(u(t,),x, (tz) 


= Uph(Uy.x, (tz) ) 


| A 


U yh (uy »M) 


< 0 because h(u,,M) < h(L(M),M) = 0 
contradiction: 


Hence u(t) < Uy: 
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Similarly we can show that if u, < L(M) then 


0 


ua b(M)™ for vali” t aE 


We also observe that no trajectory of (4.1) leaves the octant R° 


since existence and uniqueness of solutions is assumed to hold 


because if u 0 then wu = 0 


tad 
ul 


0 then Xt = 0) 


! = 
7 0 then x5 0. 
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Results mentioned above can be put in the following form. 


Treoren<4..1 #0’ bf the fassumptironsi® Gy=s()i\tand i (Xrwige dressatistied 


then 


R = {(u,X,5X,): OP SF L(y) 3707 <"x 


sea region of stability for’ (4.1). 


Proof: We have established that any trajectory of (4.1) starting in 


Reystaysuin wR» hence the sesulé. 


Aco. Equilvbprrum States. 


In this section we list all the equilibrium states of the system 


CA Le 


E,: (0,0,0) is always an equilibrium state. 


The, conditions ..(ii),, Giv)) |) jandaye(4) aegave the following equilibrium 


states 
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Depending upon the number of intersections of the curves 
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L(x,) and 
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we shall have various equilibrium states in the U-X) plane. However, 
we can make some further assumptions on the functions h(u,x,) and 


8, (x,.W), Sovas LO guarantee the existence of a unique equiliprium 


state, interior to the U-X) plane, 4.e: 


ES: (5x6 0)., where u, 5 are such that 


L(K, (u)) =u 


| 
x § 


and Kk (ECe)) = L 
Prom ©? (i) "did? °(Gv)ih we. sgetitthe equilibrium state 


Eg: €i3(0), 0°.K., ) 2 
Again depending upon the number of intersections of the curves 
8, (x, 0) = 4, (x, >x,,9) 


and 85 (x,) = QV (x >X5) 


one can have several equilibirum states. Let us assume that X19X5 


solve the above equations, then the corresponding equilibrium state is 
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EO: (0,x. 5x5) 


/ 


Finally,,the equilibrium state,of greatest interest is the one 
wageh lies interior to the first octant {(u>0, x > 0, X54 > 0). Any 
equilibrium state of this kind will be obtained by solving the 


following system of algebraic equations: 


I 
oO 


h(u,x,) = 
8, (x, .u) = ay Xy>x,-4) 


85 (x5) = 4, (%1>%X>) 


Foom "=(2)" “we"get ~u'> L(x)), so that the existence of such equilibrium 
depends upon the solution of simultaneous equations 
aU trad cul) wet cote Boppy: BO B60) 


8, (x5) = 45 (%,>X5) 


Xx solve these equations, then we have the 


Let us assume that Xp 5 


required interior equilibrium state 


E.: (u,xX,.X,) where u = L(x,). 
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M(u,x,5X,) = 
h(u,x,)+uh,,(u,x, ) pen oe 0 
[aig, (x, ,u)-q, (x,,x,,u)} 
0x1 187), (XU) -4),, (xy .X5,u)} pp ee Saw bie WAST Is 


Tie SL 


8, (x5) -4, (x, >X5) 
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Now we consider various equilibrium states separately. 


E,(0,0,0): From. (4.2), we can compute the characteristic 


equation for Ey: As mentioned in Chapter II, it is given by 
h(0,0)-A 0 0 
0 ag, (0,0)-A 0 = 0 
0 0)-A 
0 g., (0) 


since, “h. (0,0), g, (0,0), g, (0) > 0, this equilibrium state is 
unstable in all directions. So that all populations grow near Bo 


This also shows that all the populations cannot go to zero, simultaneously. 
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E5058, (0) 00): The gcharacteristic-equation tom 9£..) 15 


h(O,K, (0))-a 0 0 
ak, (0){g,,,(K, (0) ,0) aK, (0){g,, (kK, (0) ,0) AO TEMES Cuew, 
-4y4,(K (0),0,0)} —-ay, (K, (0) ,0,0)}- ; 
1 


The eigenvalues are h(O,K, (0)), aK, CD48), HK 0) ,0)-a,,, (Ky (0) ,0,0)3 


and {g,(0)-4,(K, (0) ,0)}}. From (i), (ii) and (ix) we find that 
E., is unstable in the u-direction, stable in the x, -direction and 
the stability in the x,-direction depends upon the sign of the 


expression te, (0) -q(K, (0), 0). However E is unstable. 
eae) The characteristic equation for E. 1S 


h(0,0)-A 0 0 


0 a{ g, (0,0)-4, (0,K,,0)}-A 0 Sf 
Y Se Fase (OBK ) Kot Soke Codd (0.1 )3- 


The eigenvalues are h(0,0), afg,(0,0)-q,(0,K,,0)3 and 


RD) Goseuen, Wezienay K \} (thus }near E., u-population grows, 


X,-population declines and the behaviour of xy depends upon the sign 


of the expression {g, (0,0)-q, (0,K,,0)}. Certainly E. is unstable. 
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E, (L(0) ,0,0): The characteristic equation at Ey is given by 


L(0)h, (L(0) ,0)-A L(O)h, GL (Gie.0) 0 
1 
0 ag, (0,L(0))-r 0 = 0 
0 0 8, (O0)-A 


This shows that Ey is stable in u-direction but unstable in 


directions. At low population levels of xy and X55 


X1> 
the effect o 


the mutualist is negligible. 


Next we consider the equilibria lying in the planes. 
Be ee 0) The characteristic equation at E is given by 


uh, (u,x,)-A ee 0 
ax, {1 (x, W)-4),, 0x) 0,u)) Sax re gy per vesbay ht a 
ag To ae 


0 0 {g.(0)-q,(x,,0) 3-2 
The eigenvalues are given by the roots of the equation 
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eee + =u u,Xx Xx rat x u 
r uh (u,x,) + ay Steer Diggs er tanalids 2:0 


AA, = aux. [h (u,x a) Reg 0a0 
d erie Uc vae hes ae 


“hx, (u,x,) fg) (x, ,u)-4,, (x, ,0,u) }] 


If we assume oes aoOr then AA, > OF somtnat Gn thisscase Ee is 


asymptotically stable in the u, xy plane. The stability in the 
x,-direction depends upon the sign of (g,(0)-4, (x, .0)}. In case we 
assume that ga 0, then after Rescigno and Richardson by making 


further assumptions 
uh, (u,x,) + xh. (x, ,u) See 20 


ug, (Xp) + aay Oe ec <0 


We can still have the same conclusion as above. 


E,(L(0) ,0,K 3 The variational matrix at Eee assumes the form 


L(0)h, (L(0) ,0) L(O)hy, (0) ,0) 0 
M(E,) = 0 a{g,(0,L(0))-4,(0,K,,L(0))} 0 


The eigenvalues of this matrix are 
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L(O)h (L(0) 0), afg,(0,L(0))-4,(0,K,,L(0))} and Kyl Box. (Bq) Ay, (OK) 


This shows that near E u, X, populations are stable whereas the 


6’ y 


stability in the x,-direction depends upon the sign of the expression 


a{g,(0,L(0))-a,(0,k,,L(0))}. 


E,(0,x,,x,): As mentioned earlier there may be several equilibria 


in the Xs X5 plane... Thestharacteristic equation,ein this. case.<is 
given by 

h(0,x,)-A 0 0 

ax, {811 (x, 0) OSE am Gable ce? 
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i I, whe . é ; 
one 1, where g, (x,,0) q, (Xx, 50) > 03 850%) - 4, 0%) .x5) > 93 
Zone II, where 8, (x, 0) ~ q, (Xx, 50) <70* g, (x5) - q4(X, >X5) <1) 


and Zone II, where {8 (x1 ,0)-4; Cy x, ,0) 18, (x5) -4, (x1 .x,)} ls 


Zone iit, vitae by the curves 8, (x, >0) = q(x) .x,,0), 8, (x5) = 

qo (X1 0X5); contains the points enclosed by them in R? . The 
stability of the equilibrium depends upon the sign of A542, which is 
given above. If A.A, > O then EQ 1S stable an (xX) >X5) because 


23 


ho + hs OUD sie horz <0 then Ee is a hyperbolic point in the 
plane.» In jall these cases the u-population grows near the (x1 x5) 


plane. 


E,(u,x,.x,): Finally we consider the interior equilibrium. The 

variational matrix at Ee 1s given by 

[ uh (u,x,) y Bee, on 0 
= fax, {810% >¥) oe pieuhre farrhy (4.3) 

- 43 Xp oXqo wt “ne Cree) 

y pore Te pee) 
mal (xX, >X5)} 
2 

The eigenvalues hyprAgrAg (say) of this matrix satisfy the following 


relations: 
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ta = = uh (u,x X X..uU yu 
fens 1s 2 His) sea Nah ween a be, (PS os 
7 pC ony ez, (x, x5) }] zs eon (x, ,u) 
Z 2 = 


- Vx, KL 42x, (x,>x,)] - aux,h, se x, ){8,,, (x, >¥) 


- 454, % 2X2) I 


+ Se 2 oes see ae oe ee 2 
ae a eux XU Ce )1h8) OW) Vx, Xp I B90, Op) 
Vx Hy rXQ)I-Ay x Hy r%qr4A7¢ Oy r%p)} ¥ Ba eg )tg) 0x 4) 


yy 2% pW Bp, Cy) “Ag, FQ) H (4.4) 


We now make a change of variables in the system (4.1) as follows 


Cc 
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c 
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Naa) 


ee a a (4.5) 
eee ic a 
and introduce 
vi 
Nene go) (column vector) (4.6) 


Then the system (4.1) can be represented in the form 


x= Axe F(x) (46) 
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where the matrix 


A = M(u,x, ,x,) (4.8) 


and F(x) is a column vector given by 


F,(6.n,63 U,X, 5X5) 
F(x) = F,(E.n,6; U,X] »X>) (4.9) 


cee) A aa 
P(&.n,6; U,X1 4X5) 


where 


(E+u)h (E+u,n+x,) = uh, (u,x,)é- uh, (u,x,)n 


Fi(E.n.o3 U,X, 5X5) 


Fi(E5n,6; U,X)»X>) = a (n+x,) [gy (ntx, ,E+u)-q, (ntx, ,o+x,,6+u) | 


0X, 18,,, (x) Ww) -4, (XX, WIE + Salis alae 


grat ae en ie Oneerane 


and 


F(E.n503 U,X, »X5) = (c+x,) [g,(St+x,)-G, (n+x, ,o+x,)] = aprape an ore 


+ Xo(By. (X)-ag, (XX) he 
2 : (4.10) 


The linear approximation of system (4.1) near the equilibrium state 


ss) is given by 


xy = AX (4.11) 


Theorem 4.2. Let us assme that the system (4.1) has a unique 
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incertonmsequl librium MG. xm ,%-) « R° . Also let 


yk 


(eta > Ob ae Oe anda a a, > az. 


inven ex(t) = 0, “t "2. 8G for any to is an asymptotically stable 


solution of (4.7), where a) »a,,8, are given by (4.4). 


Proof: This is proved exactly the same way as the theorem (3.5). 


teow) Persistence in. Model (4.1). 


We first define (after Gard and Hallam, 1979) extinction and 
persistence with reference to our model (4.1). The system (4.1) will 


be said to have a solution which goes to extinction provided there is 


a trajectory co of (4.1) where 


+ 5 ; 
OSM Ga ee elaad MAA RaN Dy SALET Vic MO on colle onions lye wel URE ee fo 
C 1 2 0 Tee 2 eto 0 (4.12) 


Such that at least one of u(t), x, (t), or x, (t) approaches zero as 


t approaches infinity. When no solution of (4.1) goes to 


extinction, the system is said to be persistent. 


Theorem 4.3. Let the following conditions hold, in addition to those 


mentioned in Section 1. 
(i) min{(g, (0,0)-a, (0,K,,0)), Cg, (0,L(0))-a, (0,K,,L(0)))} > 0 
(ii) min{(g, (0-4, (K, (0) 0) ), Cg, (0)-a, (x, ,0))} > 0 
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(iv) i cate oe pea ug, ,, (x1 wu) +x 
Then the system (4.1) is persistent. 


Proof: First we show that the component u(t) can never go to zero. 
Eet.=6 >.0 be given such that <6 < L(0}. Then if Co ty is such 


eieta u(t) j<06, theg 


eucwis= u(t)h(u(t),x, (t)) ¥ u(t)h(6,x, (t)) > 0 


and hence there cannot exist any sequence th > 0 9 u(t) oa ale 
thesother (cases in which a solution to (4:1) Jean go-to 
extinction are when x, (t) Ul ase t= 3 
On x, (t) mae aS) 9 Ue 7 


or both x, (t), x, (t) go to zero, simultaneously at t 


+ 
We claim that the components x, (t), x, (t) Of "C: {ecannot, go 

to zero simultaneously. Because if they do then for sufficiently 
hargew t, u(t) 1s close to, Li0) [and the per capita growth rate sot 
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and consequently x, (t) #0 as t+, contradiction. 


Next suppose that the component x, (t) >Q as t+o while 
x, (t) and u(t) persist. If this happens then for sufficiently 
large. t,. u(t) 1s. close' to its carrying. capacity =U.) X5 Ls 


close to its carrying capacity K, and the per capita growth rate of 
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ne is given by 


1 


[g, (0,£(0))-4, (0, k,,L(0))] > 0 from (i). 


But this implies that x, (t) * 0 as t > &%, contradiction! 


Finally, suppose that the species x, goes to extinction while 


Z 


us and xy persist. then, again @or sufficiently large “t, Kael ax 


as t > © because it has been shown by Albrecht et al. (1974) that 


it 


under the additional conditions (iii) and (iv) =. is globally 
D 
asymptotically stable in R’ . SO that the per capita growth rate of 


x for sufficiently large t is given by the expression 


2 


[2 100) Ae pre) 720k per om GLa) 


But this implies that x, (t) + 0 as t/> =, which gives 
contradiction. 


Conditions (i) and (ii) also guarantee that when the 
u-population is at very low level, neither xX, Nor x, goes to 
extanccron. 

Thus we have exhausted all possible cases of extinction. This 


proves the theorem. 


Note: Expressions involved Im ithe, conditions dss) Sandee? arekfor 
different eigenvalues of the variational matrix near the equilibrium 
States. For example {g,(0,L(0))-q,(0,K,,L(0) }} is one of the 
eigenvalues of the variational matrix at (EQ), as shown earlier in 
Section 4. If this expression is negative then it can be shown that 


the real parts of all the eigenvalues at (E,) are negative. In that 
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case, by considering the linear approximation of the system (4.1) 
near (EQ), we can find a strong Liapunov function for the linearized 
system and thereby establish the existence of a region wherein each 
solution of the non-linear system that is initially in the region has 
component xy which goes to extinction. 

Simi lariy i: (g,(0)-4, (x, ,0)} < 0, then there exists a region 
in R> bounded by (u,x,) plane, wherein each Solution that is 


initially in the region has component x, which goes to extinction. 


Z 


4.6.° A Special’ Case. 


In this section we consider a special case of the model (4.1), 
which has most of its features. Let the dynamics of species growth 


be given by 


De OBX. X 
il eZ 
x7 = ax, (1 -2) = T-mu (44,12) 
ne 
X5 = sx,(1- - NX, X5> 


where K,>K, are carrying capacities of the species xy and X, 
respectively. Parameters M,%,Ly5,4,8,Y555N are all assumed to be 
positive. In the absence of the mutualist (u), the model (4.12) 
represents a two-dimensional Lotka-Volterra model with carrying 
capacities, which has been discussed by many authors. Such a 


discussion appears in Freedman (1980) and our main emphasis here is to 


reconsider all different cases dealt with therein, in the light of 
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a third species, the mutualist, and to see how the behaviour of solutions 
change in the presence of the mutualist. 

First wesllist fall equilibrium states for the model (4.12). It 
is easily seen that E,(0,0,0), E,(L,,0,0), E,(0,K, ,0), E,(0,0,K,), 
E. (LotXK, ,K, ,0) and E, (Ly >0,K,) are points of equilibria in the 
feasible octant (first octant) of three dimensional phase space. If 


the condition 
(nK, -6) (BK,-1) > 0 } Cais 


holds, then there exists an equilibrium state EAs eo) in the 


Xj 9X5 plane, where 


6K, (8K,-1) u K, (nK, -8) 


x. 2s > Xe eae rae (4.14) 
il (8nK, K.-S) Z (8nk,K, 6) 


Existence of one or many equilibrium states, interior to the positive 
octant of the phase space, depends upon the solution of simultaneous 


equations 
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if (u, x, »X,) represents a solution to this system then we can show 


that Xy> satisfies the following quadratic equation 
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MLOX) - {m&SK) +8nK, K,-6(1+mLp) }x, = SK, (1+mL,-8K,) = 0 (4.15) 


Case It le es mL = 6K 


In this case, the non-trivial solution of (4.15) is given by 


m&5K) +8nK, K,-SBK., 


i m26 
“a (8-nK, ) BK 
Nank=. eows panied GA 
1 ] m2d ; 
- K,, E: 
Then ke ae [S-nx,], (4516) 
and u = Ly + @X, 


Since Ky is the carrying capacity of the species Xy> for the 


biological realization of this equilibrium state we require that 


Thus we have shown that if 1 + mL = BK, and Ky < ~ there 


exists a unique interior equilibrium E,(u,Xx,»Xx,), where u, xX and 


x are given by (4.16). 
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GCaseri le 1"+ mL) # BK, 


In this case Xx is given by 
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= Ay 
T 2m26 ( ) 


where wu = medK, + Bnk, Kk, - 5 (1+mL)). So that depending upon the relative 


valuesof the parameters occuring in the model (4.12), we can have one, 
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We shall denote the interior equilibrium by E,(u,X,.X5), fee 
exists, where x4 is given either by the first equation of (4.16) 
Oreoy. (4.17) and Le u are always given by 


K 
D u 
ihe [S-nx, ] 


~ 
i} 


(4.18) 


We shall now consider the various cases dealt with in Freedman 
(1980), as mentioned earlier, and observe the effect of introducing the 


mutualist into the system. For stability analysis, we compute the 


A : : 5 
variational matrix M(u,X) 4x5) at (u,X) 5X5) € i 
Baie +2x 2 
1 (Lo t2x,) 
maBx.,xX 2X aBx aBx 
jee 1 % | 
MiUex= ex) = of 2) ~ (4.19) 
eerez Gea Ky 1+mu 1+mu 
0 ea 6 = ey -NX 
Me K 1 
zZ 
rant | Ou 
Case A: 2B = K,> = Ky 


In. this. case. the condition 4(4.13) ‘is satistied, so. that 
Beers x) exists. It is known that in the absence of the mutualist, 
there lies a separatrix in the X19 X5 plane such that any 
trajectory of (4.12) initiating on one side of the separatrix 
approaches E. and on the other approached Eye The equilibrium C. 


lying interior to the plane is found to be unstable. Now we prove the 


following theorem for (4.12), which highlights the influence of 
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mutualism on competing species x) and X5- 


Theorem 4.4. Let the parameters m,Ly.8,5,n 7/0 (be: such that 


1 5 
Eg Sp ea 


(b) au, > 0 such, that 


Ne mL - BK., > uy > 0 


where K; K, > 0 are the carrying capacities of the species xy and 


xX, rnrespectively: Then 


2 


(1) the model (4.12) does not admit any interior equilibrium. 


(11) for every trajectory 
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x, (t) goes to extinction, i.e. 


lim, ... x, (t) = 0 
Proof: Suppose there exists an interior equilibrium (u,X,,Xx,) E es 
- - 6 - ‘ : 
then from (4.18) Xo 2A) xy < rion x, is given by “(4217), <so 
that 
5 us Yu2ameeK, 8 (14mL )-8K,) 5 
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Le 6 
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Squaring each side and simplifying we get 
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6 1 
or K, (1+mL)-8K,) < oa - 7 (mesKy +BnK, K,-6 (1+mL))} 
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contradicoion. 


To establish 


(11), we consider a function 
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the system 
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x, (t) (4.21) 
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Bom e3( 4. 1); put p30 (orm < by + &X,; so that for any arbitrarily 


small e > 0 there exists a T zat such that 


u(t) > Ly “ie COT ate 


ing 
bet -O. << Sane then 


Perm CC) ee +e mL mer etek 4 
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wa 1+mu (t) < Tem )-me 
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: : b B 1 B 
which gives (zt ~ T+mu(t) ) ale 7 nb =) 
2 2 0) 
This further implies that 


1 gers ( 8 ) 
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Now let us integrate both sides of the equation (4.21) on [Teele 


Ne ec 
V (x (ta) x (t)) t 
1 s ek OY is Poem ae 
i V(x, (7).x my | ; ve (; Ene (e ee Oke 


< “tj (2-4): «(eo Taise)s (s)fas 
tee 6) KY 1 K, 1+mL,-Hy Z 


from .(4.22) 


which gives 


4 301 ti 


(£5. 4) “naizesps 4 0 seb ap 
Me ) 


a7 


1 
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5 VG, (ix, exp -G- Tae x, (s)ds 
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froma Cay. (4.523) 


We now claim that the above relation implies lim, , x GG) t= ONE Li not, 
then from (4.23) bin iggy (x, (ep, 25 6e)) = 0, which again implies that 
x, (t) >0O as t->©2 because x, (t) remains bounded. This completes 
the proof of the theorem (4.4). 

Theorem (4.4) shows that the introduction of a mutualist might 
cause disappearance of the effect of the separatrix,) which ‘exists in 
two dimensions, and also the mutualist-competitor may win all the time. 


Actually we can show that ES is globally asymptotically stable. 


Theorem 4.5. Let the conditions of Theorem 4.4 hold. Then the 
equilibrium state Eq (Lot#k, K, 9) is globally asymptotically stable 


in the positive octant of the phase space of the variables U,X)>X>- 


Proof: From Theorem 4.4, it follows that if Gu et) 


represents a solution to (4,12) with positive initial conditions, and 


detined 2Or-alle Cc ma > 0, then 


lim, , .X, (t) =e. 


so that the w-limit set for any such orbit will lie in the U,X, plane. 
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To establish the theorem, we need to show that Ee is globally 


asymptotically stable in the positive quadrant of the u,x 


> 1 plane. 


However, this follows directly from Theorem 3.3. 


Next we consider local stability analysis of all the equilibrium 
States in Case A. 


The eigenvalues of the variational matrix (4.19) evaluated at 


i 
Equilibrium state E,(L),0,0) is found to be attracting in the 


E, (0,0,0) are? 7,056, 0 SOutLnat.. FE is unstable in all directions. 


u-direction and repelling in the other two directions. The 

characteristic values at E,(0,K, ,0) are Y,-a,6-nK, , so that near 
E., the u-population is increasing but the x1 and X, populations 
aresdecreasing:.» The characteristic roots of the matrix (4-719) “at 


E,(0,0,K.) are y; a(1-B6K,), -6, hence near Eye u-population is 
increasing but x) and X,-populations are decreasing. Considering 


E,(L)+2k K, 9); we find that ails the charactermstic roots ‘(ile. -y, 


ih 


-a, (6-nK,)) are negative, hence each population is decreasing near 


E.. The eigenvalues of the variational matrix at E, (Ly ,0,K,) are 
BK 


-Y, a(l - ae) and © =O. 8) SOntia tree Ls) stable= in, uw) eand X5 
0 


6 


directions but the stability in the xy direction depends upon the 


magnitude of parameters m and Lo: If the product mL is very 


small then Ee is locally asymptotically stable but for larger values 


of mo» Ee becomes unstable in the Xx) direction. 
r 


state Een) in the X41 9X5 plane has eigenvalues YsAy5 ; 
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ax 8x 
. i aie 
not | Ke eats | 


| (4.24) 
eae bs ea 
vox K (5-8nK, K,)x)x, < 0 


Thus -~ continues to be an unstable equilibrium. Now we consider 


equilibria interior to the positive octant. We have shown in 


Theorem 4.4 that if the product mL» is such that BK, Syloh mL then 


Enere: 1S no interior equilibrium.* Also, since S Ky it makes 


x1 = K, in the Case I, so that there is no biologically feasible 


equilibirum state*for~=1 + mL = BK,. In case BK. Soe tie ch mLo » we 


can show that there exists an equilibrium, E,(u,X,,X5). The 


Variational matrix at any such equilibrium 1s given by 


[ -y ae 0 
iG ioe Nacyat tase, RemRell 
1 eis as (ana ie K) 1+mt 
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4 2 

0 = HX - — 
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which can be written as 


where 
. 
a = 
1 
a 4 
2 
and ae a= 
3 


Rewriting a. 


YOX) X, 
K,K, (1+mu) 
YOX4 X, 
KK, (1+mu) 
OX, X, 


KK, (+mu) 


ee a ,t—i—™l 


1+mu 


St oe ey 
fy x ; 8X, xX, ‘ m&BOX, xX, 
SES 


we have 


K, (1+ma)* 


£59 2 ‘ 
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[6 (1+mL 
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(4.25) 
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Next we compute 


Ox 6x ax. 6x adx_xX aBnx. xX aBym&x_x 
-a al ye 2 a)h( lS Ay eS ee 


jen <3 K K K zs 
1 2 Lopes ie? Pe (1+ma)? | 
YOX, X, os 
—; [6(1+mu) -BnK, K, (1+mu) -m268K, x, ] (4.28) 
K_K.(1+mu) 
2 
or 
a1a, - a. = bio + boo + i. 
where 
\ ae a a5). Sr 
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As mentioned earlier Ee exists Ji ale + mL < BK, and as 
shown in Theorem 4.4, we can establish that the postive root in the 


formula (4.17) is not admissible, so that E,(u,Xx,»X5) is given by 


2 
‘ u- u?-amesx, +6 (8Ky-1-mL,) | 
a 5) 
1 2m26 (4.30) 
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1 “2 6 (1+mL) . 
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We can show that E. is an unstable equilibrium. For this let 


us consider, relation’ (4.27). for a. 
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Using the Routh-Hurwitz criterion, we find that Ee is unstable. This 


proves the following theorem. 
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GBs we v4m&6K) +5 (8K,-1-mL)) 


3 


where Kl» K, > 0 are the carrying capacities of the species x, and 


7 respectively and wu = m&SK, + 8nK, K, - 6(1+mL)). Then there 


exists a biologically feasible interior equilibrium state E, (u,X, »X5), 


X 


which is unstable. 


The table Ty: given below summarises main features of the 


Case A. 
1 6 
Th: 8 < K,; n < Ky 
equilibria u-direction x-direction 
E, (0,0,0) unstable unstable unstable 
E, (Ly ,0,0) stable unstable unstable 
E,(0,K, ,0) unstable stable stable 
E,(0,0,K,) unstable stable stable 
E. (Lo+ 2K, ,K, 9) stable stable stable 
ae araete unstable if 1+mL,>8K, 
eels pe) unstable Hyperbolic point in the X1X5 plane. 
1+mL,>8K,, there is no interior equilibrium 
BS ee eo) tegita o Ee exists and is unstable. 


The condition (4.13) is satisfied, so that there exists an 


equilibrium FeO tee in the X1 9X5 plane. It is known that in 
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the absence of the mutualist the interior equilibrium in the 


We find that in this case, the mutualist does not seem to 


cause any change in the stability behaviour. 


various equilibria is listed below in Table T 


equilibria 


E, (0,0,0) 


Eat h.,0,0) 
E00, K,,0) 
E,(0,0,K,) 

E. (L)t2K, »K, -0) 
B.(L5,0,K,) 


E_(0,x, 5x5) 


unstable 


stable 


unstable 


unstable 


stable 


stable 


unstable 


stable 


u-direction x, direction 


unstable 


unstable 


stable 


unstable 


stable 


unstable 


stable 


stable 


B° 


Local stability of 


x,-direction 


unstable 
unstable 
unstable 
stable 
unstable 
stable 
stable 


stable 


ieee 


plane is a stable attractor in the positive quadrant of the plane. 


The verification of the above list follows easily from the discussion 


of the Section A, except for the last case. 
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Now we need to show that ay = OF a. 2 0” and aja, - a. > 0, where 
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Thus we have shown that a, > OR a. >. On*-and a,a, - az > OF ina s 
justifies the last line of the Table T. 

This analysis shows that the introduction of the mutualist into 
the competitive system does not affect the stability behaviour of the 
interior equilibrium states; it only changes the equilibrium populations, 


Blologtcally, it means that tne errect Orethepspecres xX, on the species 


x is already so small that lowering it down further by the presence 


i 
of the mutualist does not create any drastic change in the behaviour 


of the system. 


6 
Case Gs < K Kgs fF ; 
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Inethisycasezthe condition~ (4715), Jis not satisfied, so that 


5 does not exist, i.e. there is no equilibrium state interior to the 


positive quadrant of x X, plane. 


ee 
We can show, as in the previous cases, that in this case 

E, (0,0,0) is unstable in all directions, E,(Ly,0,0) is stable in the 

u direction but unstable in other directions, E,(0,K, ,0) is stable in 


the x1 direction but unstable in the other two directions, E,(0,0,K,) 


is unstable in the u-direction but stable in the x and x directions, 


1 2 
E. (Lot KK, ,0) is stable in the u and x4 directions but unstable 
in the X, direction, E, (Ly 50,K,) is stable in the u_ and Xo 
directions but the stability in the x, -direction depends on thessign ot 
BK, 


alc? : ES does not exist. Here again, we see that the mutualist 


will affect the stability of Ee: 

1fnak = mL < BK,, Ee is stable in all three directions and if 
Ss mL > BK, Ee becomes unstable in the Xx) direction. Also as 
mentioned earlier, the competitive system alone does not admit any 
interior equilibria but we shall show that by the help of the mutualist 


it is possible to have interior equilibria. The interior equilibria 


are given by (4.16) and @(4-2m 


Case C,: 1. + mL, > BK,. 
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where u = m&dK, + BnKk, kK, - 6 (1+mL,) 
K 
= ae = 
Fi ty oa Pi eae 
u = Ly + UX, 
It can be shown that Xx < K) < = By a similar calculation, as done 


in Case B, we can show that EQ» in this case, is asymptotically stable 


(locally). 


Caserc: ls + mi. < -6K 


Z 0 2" 
Prom (4.17) 
ut Ne EN ek tenes 
et 2m26 


To make this feasible we impose the condition 


u > f4m2dK, +6 (8K, -1-mL5) (4.32) 


Now, we are going to show that the equilibrium corresponding to 


ut yur -4m2sK, +6 (8K,-1-mL)) 
c ee ee eee @..33) 


3 2m2é 
is always asymptotically stable. For this we shall establish that ay, 


a and a a, - az are positive, where 21 2a5,42 are given by (4.26). 
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From *(4.27,) 


x1 Xp Z 2-2 
a = oS Se ee 6 se a ae a a hy ps =, ° 


ae? 
KK, (1+mu) 


yax.X 
12 Dao cy eae ars . 
fem? x} + (L+mLy){ misk oyu -4m2 5K) 9 (8K, -1-mL, ) 


= = 7 1 
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Wotmnee (4-55) and) (4217). 
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YOX, 5X, 
bore Wise 8 
2 2 
We notice that if DL > 0 then b, > 0, which makes aja, - az an 
sO, that let us. consider by again. 
x PX ae : X1X, - 
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*m£8K, K,nx, ] 2 0-9 pausing (4.33) “and (4434). (43a) 
Next consider 
5 Three = x if 
(1+mu) Bnk, K, 6 (1+mL,+m2x, J Bn KK, 


= m26x_ + 6(1+mL)) - 8nKK 


1 Z 


m&6K,+8nK,K,-6 (1+mL,) 
" $ ee 0 1 2 
5 + om ke -4m2oK, 6 (8K,-1-mL)) 


+ 6 (1+mL,)-8nK, K, from. (4.53) 


m&SK)-8nKk,K,+6 (1+mLp ) 


fied 2 ee es 1 i?- 
5 +> 4m2oK, 6(8K,-1-mL,) 


(4.38) 


Since we are interested in equilibria for which Xy < K)> so that 


which implies that 


= <*K or mesk, + gnk,K, - 6(1+mL,) < 2me6K 


2m26 1 1 2 y he 
m&dK, - BnK, kK, + 6 (1+mL,) > 0. Hence from (4.38) we find that 
6(l+mu) - BnK,K, > 0. Now combining (4.37), (4.38) and (4.36), we 


show that by > 0. As discussed earlier Dy >0> b, >» Os «and be 
is always positive, so that using the Routh-Hurwitz criterion, the 
equilibrium is asymptotically stable. 


Next we consider the equilibrium 


cc) ere 
uu -4m26K, +6 (8K,-1-mL,) 


e K, - 
Se ie eee 2m26 » X= Ze [6-nx,]- 


Z 
Similar computation, as done in previous cases, shows that for this 
Hence the equilibrium state is unstable. 


equilibrium az < 0. 


A summary of the Case C is given in Table TO 
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equilibria u-direction x-direction 


x, -direction 


E, (9,0,0) unstable unstable unstable 
E,(Ly,0,0) Stable unstable unstable 
E,(0,K,,0) unstable stable unstable 
E,(0,0,K ) unstable stable stable 
E. (Lot&K, ,K, 0) stable stable unstable 
Perales dat at 1+mL)<8K., stable eae 
SxaGe v7 2 : iG 1+mL)>8K., unstable 
E7(0,x, x5) does not exist 

. : 1+mL,?8K, asymptotically stable. 
Eg(u,x X5) 1+mL)<8K., one stable and another unstable. 
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In two dimensions (i.e. in the absence of the mutualist) it 


as known. that all solutions initiating in: thesintertor of the first 
quadrant tend towards the equilibrium on the X,-axis se’. (K, ,0). 

We shall show that in this case, there is no interior equilibrium 
in the positive octant and also that all solutions initiating in the 
interior of the positive octant approach E, (Lot+kK, ,K, 0). We state 


the following result. 


1 
Theorem 4.7. Let the parameters 8,65,n > 0 be such that K, < a 


and aa K)> where K,>K, > 0 are the carrying capacities of the 
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species x) and X, respectively. Then for every trajectory 
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x, (t) goes to extinction and the equilibrium state E. (Lo* KK, 0) 
is globally asymptotically stable in the positive octant of the phase 


Space oF the varitablies*” wu 7X. 5x 


ae? 
Proof: The proof runs parallel to the Theorems 4.4 and 4.5. We can 


show that there is no interior equilibrium in this case, as proved 


in Theorem 4.4. Then by considering a function 


: -l/a 1/6, 
V(X) 5%) = Oy) (xX,) 5 3 Xp oX_ > 0, a > 0 


and computing its derivative along the solutions of (4.12), we prove 
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i = 1 bounded. 
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The global asymptotic stability of E. is established exactly 
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the same way was in Theorem 4.5. This completes the proof of the 


Theorem 4.7. 


In this,case,theirole.of the jmutualistiisnot,.that important. 
Conditions are such that the species xy is already in a dominating 


position, so that the mutualist just enhances this dominance. Other 


equilibria and their stability as listed in the Table T 
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E, (0,0,0) unstable unstable unstable 
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a7, berrodic7solutions . 

In this section, we use the Hopf-Bifurcation theorem to 
establish that an ecological system, which is modelled by the general 
system (4.1), can exhibit small amplitude oscillations if the 
functions N,8) 58504) and q, satisfy a certain set of conditions. 

Let us assume that (U,X, 5X5) is an equilibrium state for the 


system (4.1), interior to the positive octant of the three 


dimensional phase space. u, Xj» X, are given by 
u = L(x,) 
g, L(x) = aCe, 355 EGaDe (4.39) 


g,(x,) = 4,(x,,x,) 


The characteristic values of the variational matrix, evaluated at 


(u,X, »X,) are given by the following equation, which is obtained from 


(4.2) 
uh, (u,x,)-A en 0 
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Big eae S Tix Xerxes 
Aaa ea 
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0 ME q5u 0 (Xz, X5) hee Gre 
Dae x. se qx (x, >x,)} r 
when expanded this assumes the form 
ie - aja? = a, = az = 0 (4.40) 
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We shall now prove the following Hopf-Bifurcation Theorem. 


Theorem 4.8. Let (u,x,.x,) be an interior equilibrium state of the 
system (4.1), lying in the positive octant of the phase space of the 


variables U,X) 5X5. Also, let the following conditions hold 
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parameter) passes through > there appear small amplitude periodic 
solutions of the system (4.1), bifurcating from the equilibrium 


state (U,X, 4X5). 


Proof: First we observe that at a = Oy » the characteristic equation 


has one real and a pair of imaginary roots because at a= a 


20, 7a. > Oe Eromim (i): © and 


r a i re 
(4192783 aca, = bj + 55% + bz = O,etromeae( tig and: (4.41). 


Next we shall show that the eigenvalues cross the imaginary axis with 


non-zero speed at a = Oy: Toido\this, let.us assume that for values 


of a near the eigenvalues are of the form 


Og» 


where Ap rAgors are real numbers. In view of this, the characteristic 


equation will look like 
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Now comparing equations (4.40) and (4.45), we get 


a, = -(A3+2))) 

,MOBe, aFureegigy § (4.46) 
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2 
(a,a,-a,) + 2h {a,+(a,+2r,) at (4.47) 


Since we want to know how rt varies with respect to the parameter 


a, we differentiate relation (4.47) with respect to 


: dn, , 
da 64143783) + 2 Gg tajt(a,t2a,)°} 


da, da, di, 
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We know that at a = Oa» Casa, = 0, hence 
dns Z 
(2b, a)+b.) + \eaaa ta,tai} iy = 0 
A=a 0 
0 
da 2b. a+b 5-4 b 
1 v0 72 Peo les : 
ae da / a=a a 2 Es 2 # 
ES) a etSiigee=a, 


This shows that as a passes through Qy » the eigenvalues corss the 
imaginary axis transversally, ie. with. non-zero: speed. 


Also, from (4.46) we observe that 


Now application of the Hopf-bifurcation theorem, mentioned in 


Chapter II, proves the theorem. 


Note: Critical value or of the bifurcation parameter a is a root 


of the following quadratic equation 
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where b,,b,,b are given by (4.42), (4.43), (4.44) respectively. 
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Le by = 0, then we have a unique a 9 1 Je 
Jhgos ks 
0 b, 


in which case we will require b,*b.z <0 Sl0#eget a feasible bifurcation 
Walue. But if by # 0, then depending upon the relative magnitudes 


ana Sagns: Of D..,b:,b we will have one, two or none of such values. 


AO ges 


In general 
-b,+ yb5-4b >. 
a caked A ip ia (4.47) 


Since G iS taken ito be positive, only positive roots in. (4:47) ‘are 


admissible. 


f782 4 opecial Case, 

Here we consider a specific model, which incorporates a cost 
associated with the help to the mutualist by the mutualist-competitor. 
The model considered in the Section 6 does not have this feature. 

This would be applicable to those biological systems in which the 
growth rate of the mutualist-competitor is hindered by the presence 

of the mutualist, when the other competitor is not around. Help to the 
mutualist-competitor by the mutualist comes only if there is another 


species which competes with the mutualist-competitor. We consider the 
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the following model 
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Bere aaraneters a,8,y,5,n,a,0,m,€,L,K, ,K, are all positive. 

In this section, we shall find specific conditions, in terms of 
the above parameters, for the existence of periodic solutions of the 
system (4.48), using our Theorem 4.8. For this we would be interested 
in interior equilibria only. 


The interior equilibrium, if it exists, is obtained by solving 


the equations 
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Thus we can expect up to three equilibrium states in the interior to 
the positive octant. The existence of any such equilibria depends 
upon the relative magnitudes of all the parameters. 

Now let us assume that there exists at least one interior 


equilibrium E*(u,X,.X5); where x, satisfies the equation (4.49). 


1 


onde, 'x. ‘are: given sby 


1s 


i= L(1+6x,) 
K (4.50) 
A 2 - 
xX = a ee 
We further assume that Xy < ~ (4.51) 
so that x, > 402 


Comparing model (4.48) with the general model we find that 
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Now we find various partial derivatives, required to compute the 


stability conditions. 
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uf 
h (ux) oe ie E ’ h (u,x,) ="ye 
g Eu 1 5 
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The characteristic equation at (u,X, 5X5), then can be obtained by 


using equation (4.40). It is given by 
hoo Aa he eek Eas a =O (4.54) 


where 2,,45,a, are to be computed from (4.40), (4.50) and (4.53). 


We can show that @,,a5,a, assume the following form in this case 
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The following result is a Corollary to Theorem 4.8. 


Corollary 4.9. Let a,24,,b),b,,b. be as given by (4.55), (4.57)- 
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2b 5% + b, # 0, then the parameter a acts as a bifurcation parameter 
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We now give numerical examples to illustrate the validity and 


practical realization of our Theorem 4.8. Consider 
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Garciations. at. E> is a: xh show that 
ay = 00962 +". 05490_> 0 
az = .0008a > 0 
by = -.0002, b, = -.0006, b. = .0897 


and the bifurcation value of the parameter a comes out to be 


approximately 19.73. 


Example 2. We construct this example by using perturbation techniques, 


First we find conditions which guarantee the existence of a periodic 
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Orba tein the vu ~xX, plane, Then we introduce the species x, and 
choose parameters in such a way that even a little above the ux, 


plane, growth rates of species U,X) »X5 remain close to zero 
and the periodic orbits of the plane now become three dimensional 
perrodic orbits. 


In two dimensions system (4.48) becomes 
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E*(u,x,) where U = 6, xX. = ),45 an equilibrium state for | (4.02) 


iM 
and further that we get periodic orbits (Hopf-bifurcation) whenever 
G@e= toy 2.0: 


We now consider the three dimensional model (4.48) and choose 


the remaining parameters as 


Kea. O20. 27308 noz, 0.2639, m= 1 and Beas 501. 


We can show that at u=6, x, =l, X, = 0.0365, growth rates of 
species U,X, 4X, are OCLOre Hence correct up to four decimal 
places £*(6,1,.0365) is an equilibrium state for (4.48). 
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so that there are two values of a, i.e. .0308 and 16.1565 at 
which aja,-az = O and a, > 0, a, > 0. Thus all the hypotheses 
of the Theorem 4.8 are satisfied and we shall have perturbed periodic 


solutions in three dimensions. 


4.9. Summary. 

In this chapter, a competitor-competitor-mutualist system has been 
modelled and analyzed. Conditions for the existence of equilibria were 
given, and the stability of these equilibria determined. As well, 
conditions for the existence of periodic solutions have been determined. 

The discussion of Section 4.6 shows that the mutualist will play 
a very important role in an ecosystem modelled by model (4.12). For 
example in Cases A and C, in which the inhibitory effect of the species 


x on the species Xi is high (i.e. 8 large), the mutualist reduces 


2 
the effectiveness of the competition coefficient 8 and thereby could 
cause the reversal of stability of the equilibrium state Ee in both 
the cases and reverse competitive outcome, provided the parameter m 
is sufficiently large. Also, we know that in the absence of the 
mutualist, the competitive subcommunity of our model (4.12) does 

not admit any interior equilibrium in Case C, but the introduction 

of a mutualist into the system always gaurantees an equilibrium 

state. This shows that the mutualist could change competitive 
exclusion to. coexistence. 


In Case D, where nonequilibrium conditions hold, the effect of 


the mutualist is to cause the competitor X, to-go torextinction. 
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In this case, the mutualist u and the competitor x, approach 


1 
equilibrium conditions. 

In Section 4.8, by considering a model, which incorporates a 
Cost, CosthemnuCualist-competitor of providing direct benefits to the 
mutualist (model of Section 4.6 does not have this feature), we have 
established the possibility for the existence of periodic fluctuations 
in populations of the community. Here, by assigning relative 


numerical values to the parameters of the system, we have demonstrated 


the feasibility of such oscillations. 
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APPENDIX 


In this appendix we shall perform the stability analysis of the 
periodic orbits, guaranteed by Theorem 3.9 of Chapter III, using the 


centre manifold theory. We rewrite the system (3.39) 
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Then we make a change of variables given by 
u- u* = Uy 
Kak y= xy (2) 
et <a a 
where (u*,x*,y*) is the interior equilibrium state for the system 
(1) and is given by (3.40). In terms of the new variables, system 
(1) becomes 
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where P(x ¥z)> PF, (uu) .X) +1) and P.(u, 5X, 5¥,) are same as 
given by (3.48) 

According to Theorem 3.9, the zero solution of (3), will 
bifurcate into periodic orbits when the bifurcation parameter a 
passes through a critical value Oo - For stability analysis of these 


periodic orbits, the_first= Step is to transform the variational 
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To effect the canonical transformation of the matrix A, we need to 
find out eigenvectors corresponding to the above eigenvalues. We 


shall first find an eigenvector corresponding to the real eigenvalue 


2. If we denote this by the column vector B_ then 
AB = 1B (7) 
Peo = [Bi eB Tt then from (4) and (7), we get the following 


equations to determine BL BAB. 
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Similarly we can find eigenvectors corresponding to the pair of 


imaginary eigenvalues in (6). For this we have to deal with the 


following augmented matrix 
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Thus the eigenspace of the eigenvalues of the linearization matrix A 
is spanned by vectors given by (9) and (10). Now we consider the 


following matrix T, whose columns are the above mentioned eigenvectors 
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The matrix T has the property that a new set of variables VioVoeVs 


(say) and the old variables Uy >X15Vyz> connected by the relation 
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are such that the system (3) transforms to the required canonical 
form in the variables Vy oVo0Vz5 which are referred to as canonical 
variables. The matrix of the transformation (i.e.T) is a non- 


singular matrix, which could be checked by evaluating the determinant 
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We can show by expanding these expressions that Gi (Vy 5V55Vz); 
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Observing the symmetry in (31) and (32) we find that the 


corresponding partial derivatives of ¢ (2) can be obtained by making 


the following changes 
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Then the periodic orbits occur for a> a, and are asymptotically 
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stable. But if the above expression is positive then the periodic 


Orbiles Occur fOr ‘a. < Oo and are unstable. 
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